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Recent lattice data have reported an infrared suppressed, positivity violating gluon propagator which is nonva- 
nishing at zero momentum and a ghost propagator which is no longer enhanced. This paper discusses how to 
obtain analytical results which are in qualitative agreement with these lattice data within the Gribov-Zwanziger 
framework. This framework allows one to take into account effects related to the existence of gauge copies, by 
restricting the domain of integration in the path integral to the Gribov region. We elaborate to great extent on a 
previous short paper by presenting additional results, also confirmed by the numerical simulations. A detailed 
discussion on the soft breaking of the BRST symmetry arising in the Gribov-Zwanziger approach is provided. 



I. INTRODUCTION 



As is well known, quantum chromodynamics (QCD) is confining at low energy. Confinement means that it is impossible to 
detect free quarks and gluons in the low momentum region as quarks form colorless bound states like baryons and mesons. Even 
if one omits the quarks, pure SU(N) Yang-Mills gauge theory remains confining as gluons form bound states known as glueballs. 
Hitherto, confinement is still poorly understood. There is a widespread belief that the infrared behavior of the gluon and ghost 
propagator is deeply related to the issue of confinement and, therefore, these propagators have been widely investigated. In this 
paper, we shall use the following conventions for the gluon and the ghost propagator, 

(a;(- p )A*(p))=8"W)(§,v-^) , (c a (-p)c»(p))=8 ab g(p 2 ). (1) 

Until recently, lattice results have shown an infrared suppressed, positivity violating gluon propagator which seemed to tend 
towards zero for zero momentum, i.e. £> (0) = 0, and a ghost propagator which was believed to be enhanced in the infrared I UBl, 
Q (k 2 « 0) ~ 1 /k 2+K with K > 0. Different analytical approaches were in agreement with these results (e.g. I3ll4ll^l^l^li^l9 JlOir ) 
to quote only a few). For instance, several works based on the Schwinger-Dyson or Exact Renormalization Group equations 
reported an infrared enhanced ghost propagator and an infrared suppressed, vanishing gluon propagator, obeying a power law 
behavior characterized by a unique infrared exponent, as stated bya sum rule discussed in JH 3,120. The infrared propagators 
have also been studied from a thermodynamical viewpoint in [11]. Also the Gribov-Zwanziger action predicts an infrared 
enhanced ghost propagator and a zero-momentum vanishing gluon propagator ||8t |9fl. This action was constructed in order to 
analytically implement the restriction to the Gribov region Q., defined as the set of field configurations fulfilling the Landau 
gauge condition and for which the Faddeev-Popov operator, 

M ab = -d,(d,8 ab +gf cb A^ , (2) 

is strictly positive, namely 

Q. = {A% 9X = 0, M ab >0}. (3) 



'Electronic address: david.dudal@ugent.be 
^Electronic address: gracey@liv.ac.uk 

^Electronic address: sorella@uerj.br. Work supported by FAPERJ, Fundacao de Amparo a Pesquisa do Estado do Rio de Janeiro, under the program Cientista 
doNossoEstado,E-26/l00.6l5/2007. 
^Electronic address: nele.vandersickel@ugent.be 
"Electronic address: hemi.verschelde@ugent.be 



2 



The boundary, 3£2, of the region £2 is called the (first) Gribov horizon. This restriction is necessary to avoid the appearance 
of Gribov copies in the Landau gauge related to gauge transformations 01. However, this region £2 still contains a number of 
Gribov copies and is therefore still "larger" than the fundamental modular reg ion (FMR), which is completely free of Gribov 
copies. Unfortunately, it is unknown how to treat the FMR analytically I 1% [l3l Il4l Il5ll . 



However, more recent lattice data Ha [171 LUa, UM at larger volumes display an infrared suppressed, positivity violating gluon 
propagator, which is nonvanishing at zero momentum, , i.e. £> (0) + 0, and a ghost propagator which is no longer enhanced, 
Q (k 2 « 0) ~ l/k 2 . This implies that the previous mentioned analytical approaches are not conclusive. It is worth pointing 
out that, recently, the authors of j20l l2lll have obtained a solution of the Schwinger-Dyson equations which is in agreement 
with the latest lattice data. Furthermore, as we have shown in a previous work 12211 . this agreement can also be found within 
the Gribov-Zwanziger approach. In this framework, we have added a novel mass term to the original Gribov-Zwanziger 
action. This new term corresponds to the introduction of a dimension 2 operator. We recall that by including condensates, 
which are the vacuum expectation value of certain local operators, one can take into account nonperturbative effects which 
play an important role in the infrared region. During the course of the current work, it shall become clear that we also 
have to add an additional vacuum term to the action, which will allow us to stay within the Gribov region £2. The previous 
paper [22] only gave a brief account of the consequences of adding the mass operator to the original Gribov-Zwanziger ac- 
tion. For this reason, here we shall present an extensive study of the Gribov-Zwanziger action with the inclusion of the new parts. 

The purpose of this paper is fourfold, and it is organized as follows. The first aim, discussed in section II, is to give a detailed 
proof of the renormalizability of the extended action. Therefore, we first present an overview of the Gribov-Zwanziger action, 
Sgz, in the Landau gauge which implements the restriction the Gribov region £2. Next, we add the local composite operator 
S,„ = \ j d 4 xA 2 to this action and we prove the renormalizability of this extended action, Sgz + S m . Subsequently, we show 
that by adding another term, Sm = M 2 J d 4 x (qkp — coco) + 2 ( N iN ^ c,X 2 , the renormalizability is not destroyed. In summary, 

section II establishes the renormalizability of the action Sgz + S m + Sm- The second aim, investigated in section EH, is to 
demonstrate that this extra term enables us to obtain propagators which exhibit the desired behavior. In particular, the tree level 
gluon propagator is calculated explicitly and the ghost propagator is determined up to one loop. Both the ghost and the gluon 
propagator are in qualitative agreement with the latest lattice results. Up to this point, we have added this mass term by hand. 
Hence, a third aim is to obtain a dynamical value for M 2 . Section IV presents this dynamical value. An estimate for the one loop 
gluon propagator at zero momentum as well as for the ghost propagator at low momenta is given. Also, the positivity violation of 
the gluon propagator is scrutinized and compared with the available lattice data. The last aim is to highlight the BRST breaking 
of the Gribov-Zwanziger action, which is presented in detail in Section V. We already stress here that it is the restriction to the 
Gribov region £2, implemented by the Gribov-Zwanziger action, which induces the explicit breaking of the BRST symmetry. 
Further, we provide a few remarks on the Maggiore-Schaden approach to the issue of the BRST breaking B23I1 . and we revisit a 
few aspects of the Kugo-Ojima confinement criterion B2411 . We end this paper with a discussion in section VI. 

II. THE EXTENDED ACTION AND THE RENORMALIZABILITY 
A. The Gribov-Zwanziger action 



We begin with an overview of the action constructed by Zwanziger 1 10] which implements the restriction to the Gribov region 
in Euclidean Yang-Mills theories in the Landau gauge. We start from the following action, 

S h = Sym + J d 4 x (b a d M A; + c*d„Dfc b ) + y 4 J d 4 x h(x) , (4) 

with Sym the classical Yang-Mills action, 



Sym = l - f d 4 xF« v F« y , (5) 



and h(x) the so called horizon function, 



h(x) = g z f c Al{U- x ) ua f ec A e ll . (6) 
The parameter y, known as the Gribov parameter is not free and is determined by the horizon condition: 

(h(x)) = d(N 2 -l), (7) 
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where d is the number of space-time dimensions. The nonlocal horizon function can be localized through a suitable set of 
additional fields. The complete localized action reads 



S = S + S y , (8) 



with 



yib b 



So = Sym + J d 4 x (b a d fl A'' + c"d M D° 

+ J d 4 x((p-3 v (d y ^+gf bm A h v <ff) -TS«d v (^(S% + gf abm A^) -g(d v W?)f bm (Dvc) b <ff : ) , 

s y = -fg J d 4 x (f ahc Ay;+r b ^t; + - {N 2 - 1) y 2 



(9) 

The fields ^<p^ e , uff^J are a pair of complex conjugate bosonic fields, while (co^ c , C0^ e ) are anticommuting fields. Each of these 
fields has 4 (A^ 2 — l) 2 components. We can easily see that the action So displays a global U (/) symmetry, / = 4 (A^ 2 — l), 
with respect to the composite index i = (jn,c) = 1, ...,/, of the additional fields ^(p^ e ,(p° e ,ro^ c ,co^ £ ^ . Therefore, we simplify the 
notation of these fields by setting 

(cp^cp- o^,<) - (9?,9f,6>?,G>?) , (10) 

so we get 

So = S YM + J d A x (fc fl 3 / ^» + c a 3 /J (^ C ) a )+ |d 4 x((pl'a v (D v 90"-^'av(D v O30 u -§(3v«f)r fo "H^vc)''(pr) ■ (H) 

Now we shall try to translate the horizon condition (0 into a more practical version [9]. The local action S and the nonlocal 
action Sf, are related as follows, 

J &A&bdc&ce~ s >' = j dAdMcdcd(pd(pdcodcoe _ ' s . (12) 

If we take the partial derivative of both sides with respect to y 2 we obtain, 

- 2f (h) = (gf^Affi) + (gf" hc A'fi b ;) ■ (13) 
Using this last expression and assuming that y + 0, we can rewrite the horizon condition 

(gf bc A;<V h ;) + ^f bc Kt;) + rfd{N 2 - 1) = . (14) 
We know that the quantum action Y is obtained through the definition 

e~ r = f d4>e~ s , (15) 



where / d<t> stands for the integration over all the fields. It is now easy to see that 

| - ° 

is exactly equivalent with equation < TT~4T >. Therefore, equation ( fT6l ) represents the horizon condition. We remark that the condition 
( [ToT l also includes the solution y = 0. However, y = would correspond to the case in which the restriction to the Gribov region 
would not have been implemented. As such, the value y = has to be disregarded as an artefact due to the reformulation of the 
horizon condition. 



As it has been proven in [9], the Gribov-Zwanziger action S is renormalizable to all orders. In the next section, we shall give 
an overview of this renormalizability, but with the insertion of the local composite operator A^A^, to extend the action further. 
Obviously, the renormalizability of this extended action S ! also includes the renormalizability of the ordinary Gribov-Zwanziger 
action S. 
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B. Adding the local composite operator A^AjJ 

If we add the local composite operator A^A^ to ([8]) one can prove l25ll that the following action is renormalizable to all orders 

S' = S + S y + S A 2, (17) 

with 

S A 2 = Jd 4 x(^A°A;-^x 2 y (18) 

with T a new source and £ a new parameter. We now go a little bit more into the details of the renormalization of this action, as it 
will be useful later. We remark that if we prove the renormalizability of the action S', we have also proven the renormalizability 
of S = So + Sy just by putting X equal to zero. We will use the method of algebraic renormalization Il28ll , Roughly speaking, 
this means that we will embed the action S' into a larger action by adding new sources, so it will display a greater number of 
symmetries. These symmetries are important as they will imply constraints on the possible allowed counterterm. The larger the 
number of symmetries, the more limitations we will find on the counterterm. This will lead to a bigger possibility to absorb the 
counterterm into the original action, thereby proving the renormalizability. In the end, we give the sources the correct physical 
values, so we obtain the action S' again. 



We shall now implement, step by step, this method of algebraic renormalization. Firstly, we introduce two local external 
sources M" 1 , V? so we can treat f" bc A'fo^ c and / a6c A"(f>f) c as composite operators just like A 2 . Hence, we replace the term Sy by 



SL 



d 4 x (M a J (D^r + V* (DM)" + *f(N 2 - 1 )) 



(19) 



If we set the sources to their physical values in the end 



phys 



= V: 



phys 



(20) 



we obtain, as requested, the term Sy defined in ©■ 

Secondly, the algebraic renormalization procedure requires this action to be BRST invariant. Therefore, we further introduce 
three extra sources N™, and r| and replace S'y + S A 2 by 



d 4 x {-m% {pw) a - g u?f»* {D,cf # + up (p^if 



(21) 



where the BRST transformations of all the fields and sources are: 







sc" 


= b a , 




= 


MS? 


= 9f, 



1 



sc 
sb a = 

S&f = 

*9i = 



Z g f»c c b c c 



(22) 



and 



sVj? = Nf* 
sr\=x, 



sM a ; = 

sN"j = 



0. 

ST = . 



(23) 
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A a 
f 


c" 


c" 


b" 


9? 


q? 


cof 




dimension 


1 





2 


2 


1 


i 


1 


1 


ghost number 





1 


-1 











1 


-1 


gy-charge 














1 


-1 


1 


-1 



TABLE I: Quantum numbers of the fields. 





"fi 




"1* 


ym 


K 


L" 


X 


n 


dimension 


2 


2 


2 


2 


3 


4 


2 


2 


ghost number 


-1 





1 





-1 


-2 








gy-charge 


-1 


-1 


1 


1 















TABLE II: Quantum numbers of the sources. 



We recall that the BRST operator s is nilpotent, meaning that s 2 — 0. We mention again that by replacing the sources with their 
physical values in the end 



K i = N u , =0, (24) 

" Iphys f Iphys v ' 

nlphys=0, (25) 



one recovers the original terms Sy + S A 2 . 
Finally, a term S ext , 



5ext — 



d 4 x(-K^D ft c) a + ^gL a r bc c b c^ , (26) 

was added, which is needed to define the nonlinear BRST transformations of the gauge and ghost fields. K£ and L" are two new 
sources, invariant under the BRST symmetry s and with 

^i Phys =n h y S =°- < 27 > 

The enlarged action is thus given by 

E = S + S s + 5 ext , (28) 

and one easily sees that the action E is indeed BRST invariant. This action now enjoys a larger number of Ward identities 
summarized as follows: 

• For the U (f) invariance mentioned before we have 



UijL = 0. 



By means of the diagonal operator Qf = Uu, the /-valued fields and sources turn out to possess an additional quantum 
number. One can find all quantum numbers in TABLE Q] and TABLEHT1 

• The Slavnov-Taylor identity reads 

5(E) =0, (30) 

with 

f .4 ( 8E 8E 8E 8E ,« 8E -a 8E « 8 ^ ai 8 ^ ai 8 £ \ 

5(E) = + ■ (31) 

• The Landau gauge condition and the antighost equation are given by 

§: = 3 ^ • C3?) 

5E , 8E 

^ +a ^ = °- (33) 
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The ghost Ward identity is 

9 a Z = Ati> (34) 



with 



,4 / § . f abc 4 ^ . m b 8 i 7?b & i ,/6i § , TT bi § 



and 



(35) 



A"i - * / dV' k ~ £V) ■ (36) 

Notice that the term A" 1? being linear in the quantum fields A^, c", is a classical breaking. 
• The linearly broken local constraints yield 



8E 
Sep" 1 ' 





5E 8E 




= gf bc A b fi U° i , 




s«) ai 1 11 




8E 5E 










8E 

"8V* 




be 


- s f abc TJ bi - o f abc A b M" 



(38) 



(39) 



(40) 



• The exact %jj symmetry reads 

%E = 0, (41) 

with 

*> - / d4t ("a, + ""4 - ■ (42) 

When we turn to the quantum level, we can use these symmetries to characterize the most general allowed invariant countert- 
erm E £ . Following the algebraic renormalization procedure M28I1 . E e is an integrated local polynomial in the fields and sources 
with dimension bounded by four, and with vanishing ghost number and -charge. The previous Ward identities imply the 
following constraints for E c : 

• The U(f) invariance: 

UijY? = 0. (43) 

• The linearized Slavnov-Taylor identity: 

«iE 6 ' = , (44) 
with 'By. the nilpotent linearized Slavnov-Taylor operator, 

/"h 4 ( 5E 5 ^LJL J^lA ifl 5 -« 5 "A 5 w-JL^ fA<\ 

SE " 7 ^SKjS 8A« + 5A« + SL" 8c" + 8c" SL" + 8c" +((> ' 8c5? +C °'' Sep? + M " 8Uf + N " W* ) 1 ' 

and 



= . 



(46) 
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• The Landau gauge condition and the antighost equation: 

?>T. C 

0, 



5E^ 
8b a 

SL SL 



= 0. (47) 



• The ghost Ward identity: 



g a L c = Q. (48) 



The linearly broken local constraints: 



8E " - ^-^g - o, 



5oT "8l# SJ " 8KJ, 



— a — - o 

5^ M SMf 



• The exact %ij symmetry: 



%£ c = . (49) 



These constraints imply that E c does not depend on the La gran ge multiplier b", and that the antighost c" and the /-valued fields 
cp", CO", (pf , can enter only through the combinations I9 H25I1 



K 


= 




Qui 


= uf- 




yai 


— V"' - 




N ai 
/< 


= <- 




M ai 

/.I 


= < 





abcfjbiqci _ gf abeybi m ci 



(50) 

The most general counterterm fulfilling the conditions d43l - d49l contains four arbitrary parameters, an, fli, «2, «3 and reads 
r = a S YM + fll |d 4 *^^+^d^+^^ .(51) 

Once the most general counterterm has been determined, one can straightforwardly verify that it can be reabsorbed through a 
multiplicative renormalization of the fields, sources and coupling constants. We also mention the renormalization factors, useful 
for later calculations. If we set = (A°, c a , c°, b", cpf , tof, (pf , c3f) for all the fields and 4> = (^, L", M"J, Nf, Vf, U™,X, r\) 
for the sources, and if we define 

go = z gg , Co = Z|£ , 

0o=z' /2 0, <I>o=Z*<I>, (52) 



one can determine 



an 

Z ? = 1 + 
Z/ = 1 +r| {^ai - —J , 

= 1 +T|(— «3 — 2a2 + 4«i — 2an) • (53) 
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These are the only independent renormalization constants. For example, the Faddeev-Popov ghosts (c",c") and the i-valued 

* 1/2 

fields (ip". co". <p7, CO?) have a common renormalization constant, determined by the renormalization constants Z g and Z A , 

Z c =Zc=Z tf =Z w = Z a =Z w =(l -r\a ) = Z; l Z A l/2 . (54) 

The renormalization of the sources (M™ ,N% 1 ,V™ ,U[}') is also determined by the renormalization constants Z g and Z A 2 , being 
given by 

Zy2 = Zm =Zn =Zy =Zu = Z g 1//2 Z A 1//4 . (55) 



Also Z x is related to Z g and Z A 2 l25ll : 



Z x = Z g Z A 1/2 . (56) 



Finally, Zj, Z/f and Z^ are also not independent as they are given by: 

Zi = Z A 1 , Z^ = Z 6 ^ 2 , Z L — Z A 2 . (57) 

C. Adding a new mass term 

1. Extended action 

We first explain the need for the inclusion of a new dynamical effect. According to the latest lattice results, the gluon 
propagator does not seem to vanish for zero momentum. This is incompatible with the actions © and ( TTTb . which both lead to 
a vanishing gluon propagator near the origin. The tree level gluon propagator in the Gribov-Zwanziger model reads l25tl 

A°(-/>MvW) - 8*0 (/) - ^) = (V - ^) , (58) 

where we have set 

A, 4 = 2^ 2 iV7 4 . (59) 

One recognizes indeed that expression d5"8l vanishes at the origin due to the presence of Gribov parameter X. In the A 2 model 
the gluon propagator is modified in the following form, 

Ki-Mw) -* g y ) ( v - ^) - 8 °> p4 +y (v - P ~f) • m 

which reveals a further suppression near the origin and thus it still vanishes. We recall here that the fields ^(p^ e ,(p° e ,ro^ e ,co™ 

were introduced to localize the horizon function 01, which implements the restriction to the Gribov-region fi. If we take 
closer look at the action (fTTT >. we observe an Acp-coupling at the quadratic level. One can suspect that a nontrivial effect in the 
tp-sector will immediately get translated into the gluon sector. For this reason, if we try to give a mass to the cp , tp-fields without 
spoiling the renormalizability of the action, we might be able to modify the gluon propagator in the desired way. Implementing 
this idea means that we add a new term to the action ( TTTb of the form Tcp^'tpf , with J a new source. If we want to preserve the 
renormalizability we have to add the mass term in a BRST invariant way. Therefore, we consider the following extended action: 

S" = S' + S^, (61) 

= Jd 4 x(-J (9?(p? - cofffl?) + pJx) , (62) 
with p a parameter and J a dimension two source, invariant under the BRST transformation 



cl 



sJ = 



(63) 
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2. Renormalizability 

The proof of the renormalizability of this action S" can be easily done with the help of the Ward identities derived in the 
previous section. Again, we embed the action S" into a larger action, 

£' = E + %, (64) 

containing more symmetries. It is subsequently trivial to check that all Ward identities d29l)-(l42l) remain unchanged up to 
potential harmless linear breaking terms and therefore the constraints d43l>-(f49b as well as the combinations d50b are preserved. 
This implies that the counterterm E 6 '' corresponding to the action E' is now given by 

E c ' = E e + E^, 

E^ = a A Ji, (65) 

with fl4 an arbitrary parameter. This counterterm can be absorbed into the original action E', hence we have proven the renor- 
malizability of our extended action. If we define 

Jo=ZjJ, po=Z p p, (66) 

we find 

Zj=Z^= Z g Z 1 / 2 , Z p = 1 +Tl(fl 4 - y - c 2 ) ■ (67) 

As the reader might have noticed, symmetries do also not prevent a term kJ 2 to occur, with K a new parameter, but we can 
argue that K is in fact a redundant parameter, as no divergences in J 2 will occur. A term of this form is independent of the fields, 
hence it would only be necessary to get rid of the infinities in the functional energy, which we calculate by integrating the action 
over all the fields 



/ 



d^e~ s " = e- w{J] . (68) 



Seen from another perspective, we need a counterterm °< J 2 to remove possible divergences in the vacuum correlators 
M(f>(p — 0x0)^.(99 — coco) \ for x — ► y. Such new divergences are typical when a local composite operator (LCO) of dimen- 
sion 2 is added to the theory in 4D. An a priori arbitrary new coupling K is then needed to reabsorb these divergences. In general, 
it can be made a unique function of g 2 such that W (/) obeys a standard homogeneous linear renormalization group equation 
ll25ll . This is a good sign, as we do not want new independent couplings entering our action or results. A nice feature of the 
LCO under study, i.e. (tpcp — coco), is that divergences °c J 2 are in fact absent in the correlators, so there is even no need for 
the coupling K here. The argument goes as follows. The Ward identitites prohibit terms in Jy 2 from occurring. Notice that 
this is not a trivial point, as naively we expect it to occur from the dimensional point of view. It is only by making use of the 
extended action and its larger symmetry content that we can exclude a term °c Jy 2 from the game. Hence, we can set J 2 = to 
find the vacuum divergence structure « / 2 , as we will employ as usual mass independent renormalization schemes like the MS 
scheme. Now, there are two ways to understand that no divergences in / will occur. Firstly, at the level of the action is easily 
recognized that the term g (3 v Co") f abm (D v c) b <$f in the action is irrelevant for the computation of the generating functional as 
the associated vertices cannot couple to anything without external CO- and c-legs. Thus forgetting about this term, the ((p,(p)- 
and (co, co) -integrations can be done exactly, and they neatly cancel due to the opposite statistics of both sets of fields. Hence, 
all /-dependence is in fact lost, and a fortiori no divergences arise. Secondly, for y 2 = 0, the action S'^ 2 _ Q is BRST invariant, 

sS'y 2 _ Q = 0. Consequently, the vacuum correlators {^((pcp — coco) ^((pcp — coco) y ^ = (^s (cpco). t (cpcp — coco) y ) = 0. Therefore, we 

have again proven that no divergences in / appear. For y 2 + 0, the BRST transformation s no longer generates a symmetry (see 

section V), hence a nonvanishing result for the correlator ^ (99 — COCo) x (99 — coco)^ or the condensate (tpcp — coco) is allowed. 

A nonvanishing VEV for our new mass operator is thus exactly allowed since the BRST is already broken by the restriction to 
the horizon. From the first viewpoint, the (9,9)- and (co, co) -integrations will no longer cancel against each other, giving room 
for /-dependent contributions in the generating functional, albeit without generating any new divergences. 

D. Modifying the effective action in order to stay within the horizon 

1. Extended action 



A very important fact is to check if it is still possible to stay within the Gribov region £2, after adding this new mass term. 
This can be investigated with the help of the ghost propagator Q (k 2 ), which can be easily read off from the Feynman diagrams 
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FIG. 1: The one loop corrected ghost propagator. 



depicted in FIG. 1, 



g ab (k z ) = 8" b g(k 2 ) = 8« b [- + - 
1 



1 1 



k 2 k 2 

4- 



N 



d 4 q (k — q) M k v 
W^\J j2K~f (k-q) 2 



V b -(l + o(k 2 )) + o(g 4 ), 



with 



a{k 2 



N g 2 f A (k-q),k v i . 
N 2 -\k 2 J (2k) 4 (k-q) 2 WW 



Going back to the original formulation of Gribov [7], being inside the region Q., is equivalent to state that 

a(k 2 ) < 1 , 

which is called the no-pole condition. In this case, the ghost propagator can be rewritten in the following form, 



g(k 2 ) 



l 



l 



k 2 l-a(k 2 ) 



4 ) 



(69) 



(70) 



(71) 



(72) 



which represents the fact that we are working at the level of the inverse propagator or equivalently, at the level of the 1PI zi -point 
functions, which are generated by the effective action F. This form is more natural, as we can now impose the gap equation dT6b . 
which is also formulated at the level of the effective action. However, in the next section, it shall become clear that the current 
action S" does not guarantee us that we are located within the region Q. as o(0) > 1. Therefore, we add a second term to the 
action, S en , given by 



$en — 2 



d(N 2 -l) 



J d d x q y 2 J 



with q a new parameter. We have introduced the particular prefactor of 2 



d(N 2 -l) 
j2g 2 N 



(73) 



for later convenience. As it is a constant term, 



is it comparable with the term (— / d 4 xA(N 2 — 1)7*) in the original Gribov-Zwanziger formulation (O. Therefore, it can be 
responsible for allowing us to stay inside the Gribov horizon by enabling a to be smaller than 1 . The explicit calculation of a 
will be done in the next section, but we can already intuitively sketch the reasoning why a will be altered. As this new term is 
independent of the fields, it will only enter the expression for the vacuum energy. However, due to the gap equation ( TT6I ). it will 
also enter in the expression of the ghost propagator (and analogously any other quantity which contains y 2 ). Recapitulating, the 
complete action now reads, 



s'" = s"+s e 



(74) 



with 5" given in equation ([6TT >. 



2. Renormalizability 

The renormalizability of S 1 " can be easily verified. Therefore, we replace S en with 

L en = J d 4 xq&J , 

with © a color singlet and BRST invariant source, s© = 0. In the end, we give © the physical value of 

d(N 2 -l) 







phys 



-Y 2 



(75) 



(76) 
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to return to the original action S'". Again, we embed the action S'" into a larger action £", 

E" = E' + Een, (77) 

with E' given by ( f64b . Firstly, as it is easily checked, the term E en can only give rise to an additional harmless classical breaking 
in the Ward identities. Therefore, all the previous Ward identities will remain valid. Secondly, we have the following additional 
Ward identity, 

8E" 

80 = (?8) 

which implies that the counterterm is independent from 0. Taking these two argument together, we can conclude that the 
counterterm will be exactly the same as before, given by d65l ). Therefore, 

snfrft = sN_ (79) 

go 8 

and consequently, no new renormalization factor is necessary, 

Z s = Z ? Z- 2 V . (80) 



3. Boundary condition 



Introducing a new parameter q, requires a second gap equation in order to determine this new parameter. We recall that, in the 
case in which M 2 = or equivalently in the original Gribov-Zwanziger formulation, we have 



o(/t 2 «0) = 1-Ck 2 , 

with C a certain positive constant, which causes the enhancement of the ghost propagator g (k 2 ) at zero momentum, 



(81) 



(82) 



Therefore, we know that at zero momentum, slowly switching off M 2 , will cause a(k 2 = 0) going to 1. It is therefore very 
natural to demand that this transition has to occur smoothly by imposing the following boundary condition, 



3a(0) 



dM 2 



0. 



M 2 =0 



(83) 



In summary, we have now two gap equations. Firstly, the gap equation tt = fixes y z as a function of M l and secondly. 



demanding that 

& dM Af =0 

discussed in sectionllVl 



will uniquely fix q. This leaves us with one free parameter, M 2 , the fixation of which shall be 



III. THE MODIFIED GLUON AND GHOST PROPAGATOR 

Now that we have constructed the action S'", by adding two additional terms 5^ and S en to the original Gribov-Zwanziger 
action, we investigate the gluon and the ghost propagator in detail. For the calculations, we have replaced the sources J and x 
with the more conventional mass notations M 2 resp. m 2 . 



A. The gluon propagator 

We shall first examine the tree level gluon propagator. In order to calculate this free gluon propagator we only need that part 
of the free action S" containing the A-fields and the <p, cp-fields. This free action reads 



d 4 x 



i(3,A«-a v A;) 2 + ^(3^) 2 +^8 2 cp; 



fg{f 



abc 



be 



f bc A aqbc ) _ M ^ae + -~ A 



m 2 .2 

V 



(84) 
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where the limit a — > is understood in order to recover the Landau gauge. The ". . ." stand for the constant terms —d(N 2 — Vff 
and 2 4 ^ = <; "f-M 2 and other terms in the ghost- and CO, co-fields irrelevant for the calculation of the gluon propagator. Next, 
we integrate out the tp- and cp-fields. As we are only interested in the gluon propagator, we simply use the equations of motion, 

dS" dS" 

r—r-r = and ^-4- = 0, which give 



We use this result to rewrite Sq, 



be = ybc \ ^2 f abc A a 



(85) 



d 4 x 



(d p A" v -d v A$ 2 - 



2a 
1 

2a 



d 2 -M 2 



gf dbc At) 



,4„2 aci 



1 



4". 

^d 2 -M 2 " 



(86) 



The last step is explained with the following relation, 



pabc fdbc ]\f$ a ^ 



(87) 



and we restrict ourselves to the color group SU (N) throughout. We continue rewriting Sq so we can easily read the gluon 
propagator 



c" 

OA — 



A ab = 



d 4 x 



-d 2 



-A a A ab A b - 

2 /i Ai /i /'V /i V 



-m 2 — 



2g 2 Nf 
d 2 -M 2 



a 



- 1 



sab 



(88) 



The gluon propagator can be determined by taking the inverse of A^J and converting it to momentum space. Doing so, we find 
the following expression 



K(p)4(-p)) = 



1 



2 2 i 2g 2 Nf 



-M 1 



sab 



p 4 + (M 2 + m 2 ) p 2 + 2g 2 AY + M 2 m 2 

■Dip 1 ) 



PvPv 



J/JV 



sub 



(89) 



From this expression we can already make two observations: 

• D (p 2 ) enjoys infrared suppression. 

• D (0) M 2 , so the gluon propagator does not vanish at the origin. Even if we set m 2 = we still find a nonvanishing gluon 
propagator, so we want to stress that this different result is clearly due to the novel mass term proportional to cptp — coco. 

In section HV Bl we shall uncover a third property, namely that £> (p 2 ) displays a positivity violation. Also this observation is in 
accordance with the latest lattice results [29]. 



B. The ghost propagator 

The observation that m 2 = does not qualitatively alter the gluon propagator, will be repeated for the ghost propagator. 
Henceforth, we set m 2 = 0, which also improves the readability of the paper. However, all calculations could in principle be 
repeated with the inclusion of the mass m 2 . 

We start with the expression for the ghost propagator. Substituting the expression of the gluon propagator, we find, 

,,2, _ N g 2 f d d q {k-q)„k v , . 
° [k > ~ W^TWJ (2%Y (k-q) 2 

2 kfjcv f d d q 1 q 2 +M 2 
8 ~W J (2%) d (k-q) 2 q 4 +M 2 q 2 +X 4 



8,uv - 



q^qv 



(90) 



13 



where we have also used equation d59l l. As we are interested in the infrared behavior of this propagator, we expand the previous 
expression for small k 2 , 

, ? x r.k u kvd-\« f d d q 1 q 2 +M 2 ,„,. 

o(k 2 *0) = Ng 2 ^—8, v / -r^-- -; h^- 2 : , 4 +0(k 2 ) 



For later use, let us rewrite o(0) as 



£ 2 {2%) d q 2 q 4 +M 2 q 2 + X 4 

d d q 1 <r+M 2 
(27t) c/ ^ 2 ^ 4 +MV + ^ 4 



7 d-\ f d d q 1 o 2 +M 2 . 9n 



1^-1 /" d <? 1 , ^2^2^-1 f d d q 1 



° (0) Ng J ' (27I)^ 4 +MV + X 4+ ^ M 7 (2n) d q 2 q 4 +M 2 q 2 +X 4 ' (92) 



Notice that the first integral in the right hand side of equation ( 1921 diverges while the second integral is UV finite in AD. 

We continue with the derivation of the gap equations as we would like to write X 2 as a function of M 2 , i.e. A. 2 (M 2 ), in 
expression ( |92l . Firstly, we calculate the horizon condition (fTol i explicitly starting from the effective action. The one loop 
effective action is obtained from the quadratic part of our action S" 

r (i) r c// 

e - r Y = / d4>e~ s o , (93) 



This time, the terms —d(N 2 — 1)Y* an d 2 ^= = q ^M 2 have to be maintained, as they will enter the horizon condition. After a 
straightforward calculation the one loop effective action in d dimensions yields, 

Setting A, 4 = 2g 2 N"f (see equation j59l), we rewrite the previous expression, 

m rl !) 2g 2 N „ 4 , , d-1 f d d q q 4 +M 2 q 2 +X 4 

= 77^^ = -A 4 + 2<;A 2 M 2 + g 2 jV— - / « in«*_ (95) 



JV 2 -1 d b 6 d J (2n) d q 
and apply the gap equation ( [TBT l. 

a^ 1 ' / M 2 , d—lf d d q 1 \ 

- 2X 2 -l+ q —+g 2 N— / ,„-,,,, = . (96) 



dX 2 \ ^X 2 6 d J ( 2 %) d q 4 + M 2 q 2 + X 4 



Secondly, we impose the boundary condition ( I83l l in order to obtain an explicit value for q. Instead of explicitly starting from 

3 

3M- 



expression d90l ) to fix q, there is a much simpler way to find the corresponding q. Therefore, we act with g|U on the gap equation 



Subsequently setting M 2 = 0, gives 

1 d-1 



^X 2 (0) d 
where we imposed d83l . Proceeding, we find 



„W^Ll 1 =0 (97) 

8 J (2n) d q 2 q 4 + X 4 (0) ' ^ ; 



d-1 , r d d q 1 



(2:t)V<7 4 + ^ 4 (0) ^ 2 (0) 
„ 9/ x 3 , /■ d 4 <? 1 1 



(2?t) 4 ? 2 q 4 + X 4 (0) 
3g 2 N 

«=ik' (98) 



which determines q at the current order. 
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With the help of the latter two gap equations d96l l and ( 1981 1, we can rephrase the correction to the self energy of the ghost. 
Combining equation 1921 and l96l l we can write 



, , d-l f d a q 1 1 M z 

o(0) = l+MYN—j ^^ ,4 +MV + nM2) ^^- (99) 

From this expression, we can make several observations. Firstly, when M 2 = 0, from the previous expression it immediately 
follows that 

a(0) = 1 , (100) 

which gives back the ordinary Gribov-Zwanziger result J3, HI US HHl ■ Indeed, from the previous expression, one derived that the 
ghost propagator, 

g(k 2 ) 



k 2 l-a(k 2 ) ' 

is enhanced and behaves like 1/k 4 , for k 2 ss 0. Secondly, when M 2 + 0, we notice that the ghost propagator is no longer enhanced 
and behaves like 1 jk 2 as already found in l22ll . which is in qualitative agreement with the latest lattice results. This behavior is 
clearly due to the novel mass term M 2 j d 4 x (cpfcpf — ©"co"). Thirdly, we see that the term in q is crucial in order to obtain a a(0) 
which is smaller than 1. Omitting this term would result in o(0) > 1 in the case that M 2 + 0. However, including this term, we 
can easily prove that o < 1 . Indeed, taking expression 1991 and replacing C, with the integral in 1981 . we find 



, 9 3 f d 4 q 1 1 M 2 N 3 , r d 4 q 1 1 



o(0) = 1+MYN-J (^4^ g4+MV + , 4(M2) -^W^ 

+ 4l 2 (M 2 ) 8 J (2n) 4 p 2 p 4 + _Mi r)p 2 + l 4X 2 (M 2 ) 8 4 8 J (2n) 4 p 2 p 4 + l 
J {2%) 4 p- 



s 2 Nl 7 ^i S7 - r — r— tt-^-tt, (101) 



3x 2 

= l ~~ 6 "] {2n) 4 p 2 (p 4 +xp 2 + l)( P 4 + l) 

with x = -J^^ > 0, hence a(0) < 1. At this point, we can really appreciate the role of the novel vacuum term (l73l . It serves as a 
stabilizing term for the horizon condition. Indeed, without the term (|73l ), we would end up outside of the Gribov region for some 
k 2 > 0, even for an infinitesimal 1 M 2 > 0. In this sense, the action S 1 " constitutes a refinement of the original Gribov-Zwanziger 
action, which is a smooth limiting case of S'". 

For later use, we can evaluate the integral in expression (l99l as it is finite. The explicit one loop value for a(0) yields 



<J(0) = 1+M 2 3g2 ^ 1 [ln(M 2 + - In ( M 2 - J - Al A] - (?£E\ 



2^A m 2 



647t 2 VM 4 -4A, 4 L V v J \ )\ V 12871/ A, 2 (M 2 ) ' v ' 

where we have substituted the value ( 1981 for q. 



In summary, we have found a ghost propagator which is no longer enhanced. So far, we have fixed A, 2 in function of M 2 and 
we have found a constant value for q. However, we have not yet fixed M 2 . This will be the task of the next section. 



IV. A DYNAMICAL VALUE FOR M 2 

Up to this point, we have only introduced the mass M 2 by hand, however it is recommendable to obtain a dynamical value 
for this parameter. We shall present two methods to find such a value. Firstly, we explain how to obtain a dynamical value for 
M 2 with the help of the effective action. However, as the calculations become too involved, we investigate a second method, the 
variational principle, and apply this to the ghost and gluon propagator, with more success. 



Notice that we must take M 2 > to avoid unwanted tachyonic instabilities. 
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A. The effective action and the gap equations 

We first explain the idea behind the method before going into detailed calculations. In the previous section we have derived the 
gluon propagator. We recall that the mass term m 2 A 2 does not qualitatively change the form of the gluon and ghost propagators, 
therefore we have put m — for our purpose. With m = 0, the tree level propagator d89b yields: 



^ = p *+^p> + 2 g 2 N f ■ (103) 

Expanding the mass M 2 as a series in g 2 , gives 

M 2 = M 2 l +g 2 M 2 +g A Ml + ... . (104) 

We only need to consider Mo, which is of order unity, as we are considering the tree level propagator. We know that at the end 
of our calculations we have to set our sources equal to zero, or 7 = M 2 = 0. If we work at lowest order, this means we have to 
set Mo = (and the gluon propagator will not display the desired behavior). However, going one order higher gives: 

Nll+g 2 M\ = 0. (105) 

The last equation might imply that Mq is no longer equal to zero, and consequently, the tree level gluon propagator will attain 
the desired form. Let us elaborate further on this aspect. 



p 2 +M 2 



1. One loop effective potential 



To implement the above-mentioned ideas, we shall first calculate the one loop energy functional. We start with the action d6"Tl ), 
whereby setting m = is equivalent with putting X = 0. We replace the mass M 2 again with the source J. In order to determine 
the one loop effective action, we first need the one loop energy functional Wo (J) which we obtain from the quadratic part of the 
action, 



From the previous expression we find for Wo (J), 



W (J) = - 



d(N 2 -l) x4 , d(N 2 -l)^ 2j , (N 2 -l) 



qX z J- 



^ 'J (2n) d 



2g 2 N g 2 N " ' 2 

We shall work in the MS scheme, and use a notational shorthand: 



(2t 



In 



P \P 



X 4 



(106) 



(107) 



/- VJ 2 -4X 4 



ml = 



J+ VJ 2 ~4X 4 



(108) 



whereby X 2 is defined in equation d59l l. Evaluating the integrals in Wq(J) gives 



This calculation is explained in detail in the appendix. 



3(N 2 -1 ) {8 
647I 2 



:^ +m 4 ln ^. +m 4 ln ^_ i/ 2 ln ^ 



(109) 



As we have determined the energy functional Wq(J), we can now calculate the one loop effective action via the Legendre 
transform of W (J) . If we define 



a(x) = 



then 



8W(J) 



a(x) = o(x) -ct c i = - 



d(N 2 -l) 2 



/ d4>((p(p — coco)e" 
Jd4>e- 5 '" 



(110) 



(111) 
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represents the expectation value of the local composite operator, — ((pep — coco). The effective action is given by 

T(o) = W(J)- Jd 4 xJ(x)a(x) , (112) 

or equivalently, as we prefer to work in the variable a, 

T(o) = W(J)- jd 4 xJ(x)(d(x)+c cl ) . (113) 

Calculating T(a) by explicitly doing the inversion is a rather cumbersome task. In most cases one can perform a Hubbard- 
Stratonovich transformation to eliminate the term J (epep — coco) from the action and introduce a new field a' which couples 
linearly to J. This greatly simplifies the calculation. However, in this case, it seems impossible to do such a transformation as 
a required term in J 2 is missing. Hence, there is no other option than to actually perform the inversion. In order to calculate 
this inversion, we shall limit ourself to constant / and 8 as we are mainly interested in the (space time) independent vacuum 
expectation value of the operator — ((pep — coco) coupled to the source J. This vacuum expectation value is given by 

~, fd4>(epcp-oxo)e^ 

a| ^° = Jd*e^ ' (114) 

where S represents the ordinary Gribov-Zwanziger action (0. As we already have calculated W(J) up to one loop is it straight- 
forward to verify that 



(115) 



d , s 13(JV 2 -1) / t ( n 1 \ 1+Vl^f 

whereby we shortened the notation by putting t = 4X 4 /J 2 . From the previous expression we find for the condensate 

3(N 2 -l) 

= --^r^. (H6) 

This is an important result, as it indicates that a nonzero value for the Gribov parameter y will result in a nonvanishing condensate 
(— ((pep — coco)) even at the perturbative level. 

We are now ready to compute the effective action up to one loop along the lines of [30]. The energy functional can be written as 
a series in the coupling constant g 2 , 

W(J) = W (J)+g 2 W l (J) + ... 

= t(s 2 ywi(j). (in) 

1=0 

As a consequence, looking at the definition dl lot , we can write 

o = So (J) +g 2 di (/) + ... 
= t(8 2 YMJ), (118) 

i=0 

where S, (7) corresponds to the ith order in g 2 (regarding J as of order unity). This is called the original series. The inverted 
series is defined as 

/ = J Q {a)+g 2 Ji{d) + ... 

= i(g 2 yjj(d), (ii9) 

7=0 

with Jj(a) the jth order coefficient. Substituting (11 19b into ( II 18b gives, 



i=0 



.7=0 



= do(Md))+g 2 (d' (jQ{d))-J l (d) + d l {Jo(d))) + ... . (120) 
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By regarding S as of the order unity and by comparing both sides of the last equation, one finds 



o = a (Jb(a)) , (121) 
_ 81 (Jo (8)) 
°o (Jo (a)) 



For the moment, as we are working at lowest order, we only need equation (1121b . We can invert this equation, so we find for 
Mo): 

Ma) = So 1 (8), (123) 

meaning that we have to solve 



a-ann I 3(iv 2 -i) r ^ jM) , I r — TTTTT j 1 V, 1+ y/i-r(Mo) ^ n94 , 

for Jo* so we can write 

7 = f(d,X) . (125) 

We immediately suspect that this inversion will not give rise to an analytical expression. Once we have found f(a,X), we 
substitute this expression into the effective action, 

T(a,X) = W(f{a,l),X)-f(d,X)d. (126) 

At this point, as we have found an expression for the one loop effective action, we can implement two equations to fix 8 and 
A,. Firstly, the minimization condition reads 

^r(8A) = o, (127) 

and secondly, the horizon condition ( [ToT l can be translated as 

J^r(8,A) = o. (128) 

We start with the first gap equation. Replacing T by equation d 1 1 31 > leads to 

9 ,^ dW dJ dJ^ dJ ^ . 

^r(o,A) = => ____o-^o cl -/ = => J = => f(a,X)=0. (129) 

Since there are only 2 explicit scales, A, and 8, present, the first gap equation can be used to express e.g. 8 in terms of A. For 
the sake of a numerical computation, we can therefore momentarily set X = 1. From FIG. |2] one can obtain an estimate 8' of 
/(8', 1) = 0, with 8' = j64tz 2 j^-j . Doing so, we find 8' w -6.28, so that 

d *^ Mx {n^y> (i3o) 

which of course corresponds to the already obtained perturbative solution dl 16b . The second gap equation ( 11281 ) must then 
consequently also give us back the perturbative solution. To check this, we first calculate the perturbative result for A- by taking 
the limit J — > in expression ( 1109t 



ro = -^^r + ^^^v-ix^). (i3i) 



Next, we take the partial derivative with respect to A- which gives, 



£ . ^(-^MfA+Mz^ti-^ | . (13 2) 

dX \ g 2 N 647I 2 V 3 V 
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fO 1 , 1 ) 




FIG. 2: A plot of f(d r , 1) in terms of S' = § 6^-^ 



The natural choice for the renormalization constant is to set /J = X to kill the logarithms. Imposing the gap equation = 
gives us, 



8 2 N = 8 
167T 2 5 ■ 

We remark that we have neglected the solution y = 0, as explained in section ULAl From 



2/-2 



and expression d 1 33b we find an estimate for X: 



1 -u a 11 N 

— , With P0 = - r T7Z2, 

poln-^ 3 l6% 

MS 



X 4 = e 44 / 15 



(133) 



(134) 



(135) 



where we have worked in units = 1. This perturbative solution is also in compliance with [25]. Now, we return to the 
effective action ( 11261 ). We first take the partial derivative with respect to X, afterwards we set = 3, we explicitly replace g 2 by 
expression ( 1134b and we use the minimizing condition ( 1130b . Numerically, we find the following value for A. 4 : 

X 4 = 1.41 , (136) 

as one can read off from FIG. [3] This is exactly the perturbative result ( 11351 ). If we calculate the vacuum energy with this value 



for X, we find from (11311 ), 



FIG. 3: The horizon function ^ for N = 3. 



£vac - Ir^e 44 / 15 . (137) 
64 Jt z 



We notice that the vacuum energy is positive. 
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2. Intermediate conclusion 



We can conclude at this point, that in the framework we have used, we recover only the perturbative solution. Unfortunately, 
at lowest order, one finds Jq = as explained in the beginning of this section, so we were unable to find a dynamical value 
for M 2 at first order. However, if we would be able to go one order higher, with Jq + g 2 J\ = 0, we might find Jq + and 
consequently the gluon propagator at tree level would attain the desired form ( 1 1 031 >- In addition, we might even discover a 
nonperturbative solution. Unfortunately, this is not as straightforward as at leading order. The main difficulty resides in the 
evaluation of two loop vacuum bubbles for the effective potential with three different mass scales. Whilst the master integrals 
are known, j3ll[32l[33ll . the main complication is that the propagator of d89l with m 2 = needs to be split into standard form 
but this introduces the masses of (1108b which are either complex or negative. In either scenario the master two loop vacuum 
bubble is known for distinct positive masses and involves several dilogarithm functions. Therefore in our case for even the 
simplest of mass choices the resulting dilogarithms will be complex as well as being a complicated function of m 2 , m 2 and A.. 
Moreover, this is prior to computing the full effective potential itself by adding all the relevant combinations of master integrals 
together. Therefore, it seems to us that whilst such a computation could be completed in principle, currently the resulting huge 
expression could not possibly lend itself to a tractable analysis similar to the relatively simple one we have carried out at one loop. 



B. Applying the variational principle on the ghost propagator and the gluon propagator 

In this section, we shall rely on variational perturbation theory in order to find a value for the hitherto arbitrary mass parameter 



2 



M 

Along the lines of J13l, we introduce a formal loop counting parameter I = 1 by replacing the action S with jS. At the same 
time, we replace all the fields <t> by v^O. Symbolically, 

S(<t>,g)^jS(yTm,g). (138) 

It is readily derived that multiplying each field with a factor of ^fl and performing an overall \ jl rescaling is the same as 
replacing the coupling g with Vig, so we can replace ( 1138b with 



S{<b,g)^S{*,y/tg). (139) 

In this fashion, the free (quadratic) part of the action is ^"-invariant, while every interaction terms contains powers 2 of \fl. The 
first order in the ^-expansion, obtained by setting I = 0, then corresponds to the free theory. More generally, the ^-expansion is 
equivalent with the loop expansion, where it is understood that we put the formal bookkeeping parameter I = 1 at the end. 

The next step is to introduce the variational parameter M 2 into the theory. This is done in a specific way: we add the quadratic 
mass term Sm = M 2 J d 4 x ((pep — coco) + 2 ^~ l ^ cjk 2 to the action, but substract it again at higher order in I, i.e. we consider the 
action 

s(^, g )^s(^,Veg)+s M -i k SM, (140) 

with k > 0. Since I = 1, we did not change the actual starting action at all. 

However, we maintain the strategy of performing an expansion in powers of £. Since the mass term is split up into 2 parts 
~ (1 — £ k )M 2 , both parts will enter the ^-expansion in a different way. At the end, we must set I = 1 again. If we could compute 
an arbitrary quantity Q exactly, the M 2 -independence would of course be apparent since the theory is not altered. However, 
at any finite order in £, a residual M 2 -dependence will enter the result for Q, due to the re-expanded powers series in I. Said 
otherwise, we have partially resummed the perturbative series for Q. by making use of the parameter I. The hope is that some 
nontrivial information, encoded by the operator coupled to 1 — £ k , will emerge in the final expression for Q,. One query remains: 
how to handle the M 2 which appears in the approximate Q,? Therefore we can rely on the lore of minimal sensitivity |[35ll : we 
know that the exact Q cannot depend on M 2 , hence it is very natural to demand that also at a finite order = 0, leading to a 



2 We recall that the perturbative expansion is one in powers of g 2 , and thus in integer powers of I. 
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dynamical optimal value for the yet free parameter M 2 . 

The described method of variationally introducing extra parameters into a quantum field theory provides us with a powerful 
tool to study nontrivial dynamical effects in an approximate fashion, yet the calculational efforts do not exceed those of 
conventional perturbation theory. 

We still have to choose a value for k. We recall that the constant term, S en , was introduced in order to stay within the horizon. 
Therefore, we want to retain this term when we are applying the variational principle. However, we are working up to first order, 
meaning that we shall expand the quantity Q_ up to first order in I and subsequently set I — 1 . Hence, taking k = 1 is not a good 
option as the constant term would vanish and have no influence. Therefore, a better option is to take e.g. k = 2, to assure the 
consistency of the variational setup with the restriction to the Gribov region. In this way, we are simply coupling the variational 
parameter M 2 directly to the theory. 



1. The ghost propagator 
We start from the expression d72| i of the ghost propagator 



a, 1 1 



and apply the variational principle on the ghost propagator near zero momentum. We have, 

o(**«0) = Ng 2d —± f^rX 4 q2+ 2 M 2 \ 4 +0(k 2 ). (142) 
V ' 6 d J (2%) d q 2 q 4 +M 2 q 2 + X 4 v ' 

As explained above, we replace g 2 — > lg 2 and M 2 — > (1 — l 2 )M 2 . Subsequently, we expand Q (k 2 ) k 2~ in powers of I corre- 
sponding to a re-ordered loop expansion. As we have calculated the ghost propagator up to one loop, we only need to expand 
the above expression to the first power of t, 

2 J- I f d d q 1 q 2 +M 2 
c(0 = Ng 2 e—— / — f-r-r , T 9 — rj ■ (143) 
W 6 d J (2%) d q 2 q 4 +M 2 q 2 + X 4 



2 d-l f d d q 1 q 2 +M 2 



As indicated earlier, setting £ = 1 gives 

«(<>) = N^-fj ^ 7 q *: M 2;; +v ■ tm 

which is exactly the same as d 142b . This expression not only depends on M 2 , but also on X 2 . However, we already know that X 2 
and M 2 are not independent variables, as they are related through the gap equation i 



M 2 7 d-l f d d q 1 

-l+c—r+g^ / — ^-j ^ — tt=0. (145) 

A, 2 d J (2%) d q 4 +M 2 q 2 + X 4 

Following the variational principle, we replace M 2 with (1 — £ 2 )M 2 and g 2 with lg 2 , expand the equation up to order £}, and set 
£ = 1 in the end. Doing so, we recover again expression ( 11451 >. At this point, it can be clearly seen that k = 1 in equation j 140b 

would cancel the effect of the constant term while k = 2 is a better choice 3 . Evaluating the integral in expression ( 11451 ), we 
find 

Ng 2 [5 „ m? , ml „ m\ , rr&\ M 2 „ A „ 



647t 2 V2 VM 4 -4A, 4 f- VM 4 -4X 4 J? J X 2 

This integral could be similarly calculated as done in the appendix, or one could start from the effective action d 1 26b and derive 
this equation with respect to A, 2 . We recall that from the boundary condition d98l l, we have already determined q = ifej. 



3 Actually, every value for k, with k > 2 is allowed. 
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FIG. 4: X 2 in function of M 2 in units Attf- = 1. 



We still require an appropriate value for fj. Therefore, we fix /r 



VM 4 - 4A 4 



which was chosen as in [25]. We 

S 2 N 



have opted for this specific renormalization scale /j which shall result in an acceptably small effective expansion parameter 
Consequently, from equation ( |134t , we find 



8 2 N 

167T 2 



llln(§ M 2 + VM 4 -4X 4 ) 



(147) 



in units of Ams = 1 ■ 



In summary, as o(0) remained the same after applying the variational principle, we can take the expression ( 1102b for o(0), 
,3g 2 N 



o(0) = l+M 2 



1 



64tc 2 



M 4 -4(A, 2 (M 2 )) 2 



In [M 2 + Jm 4 -4{\ 2 {M 2 )) z -In [M 



M 4 -4(k 2 (M 2 )) 2 
3g 2 N\ 



12871 



M 2 

) XHm 2 ) 



where A. 2 (M 2 ) is determined by the gap equation 

, , N s 2 







64tc 2 



:ln-4- 



:ln^ 



VM 4 -4A. 4 V VM 4 -4X 4 7j 2 



3g 2 NM 2 
12871 1 2 " 



(148) 



(149) 



Before continuing the analysis, let us first have a look at the gap equation. The gap equation solved for X 2 as a function of M 2 
is depicted in FIG. |U We find two emerging branches, displayed by a continuous and a dashed line. The former solution exists 
in the interval [0, 1.53], while the latter one only exists in [1.25,°o[. As the latter branch does not exist around M 2 = 0, we shall 
not consider this solution because the boundary condition d83l demands a smooth transition for the M 2 — > limit. 



We can now have a closer look at the ghost propagator or equivalently o(0). We have graphically depicted a(0) in FIG. [5] 
Firstly, from the figure, we see that a(0) is nicely smaller than 1 for all M 2 in the interval [0, 1.53]. This is a remarkable fact as 



it implies that we have managed to stay within the horizon. Secondly, we notice that the boundary condition 

indeed fulfilled, which is a nice check on our result. We can now apply the minimal sensitivity approach on the quantity tf(0). 
From FIG. [5] we immediately see that there is no extremum. However, looking at the derivative of c(0) with respect to M 2 , we 

do find a point of inflection at M 2 = 0.37 A 2 ^. Demanding = is an alternative option when no extremum is found j35ll . 

Taking this value for M 2 , we find 

a(0) = 0.93. (150) 

The effective coupling is given by 



M 2 =0 



= 0is 



8 2 N 

167C 2 



0.53 



(151) 



which is smaller than 1 . 
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2. The gluon propagator 

In order to apply the variational principle to the gluon propagator, we require its one loop correction. Given the rather 
complicated form of the propagator, obtaining the full exact expression for its one loop correction is not possible. Indeed to 
appreciate how cumbersome such an expression could be one has only to examine the M 2 = m 2 = case, 1 36], where all the 
one loop corrections to the propagators are given explicitly. However, despite this we can still achieve our main aim of studying 
the low momentum behavior of the gluon propagator corrections directly in the zero momentum limit without knowledge of 
the full correction. In ll36Tl this limit for the gluon propagator was deduced from the exact one loop computation. However, the 
resulting expression tallied with that obtained via the vacuum bubble expansion of the underlying 2-point functions. The latter is 
a much easier technique to apply and given the equivalence of the expressions it justifies its application to our case when M 2 + 0. 
Briefly one expands the 2-point functions relevant to the gluon propagator construction in powers of the external momentum p 2 . 
Though the expansion is truncated at some order such as 0((p 2 ) 2 ). The accompanying Feynman integrals are massive vacuum 
bubbles which are essentially trivial to compute at one loop. However, our situation is complicated significantly by the fact that 
there is mixing in the quadratic part of the {A^, (f^f} sector of the tree action. Therefore in addition to the gluon propagator, (70), 
we require the propagators of the remaining fields. For this derivation here we use the conventions and notation of the article 
i36Tl for an arbitrary color group, where the M 2 = problem was discussed at length. There it is evident that one has to consider 
the full {A",,(p"/'} part of the momentum space action in order to invert the quadratic sector to derive all the propagators. In the 
Landau gauge we find the set of propagators for our situation are 

Zacfabd rabe rcde^A 

(tfCdtf (-/>)> = - j^f^ + ^WTWT^ ^ (P) ' (152) 
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where the presence of 1/V2 was a key ingredient in ensuring that ghost enhancement correctly emerged in the M 2 = case, 
lf36h . Therefore we are confident that our extension here will include the previous valid analysis and therefore will provide a 
useful check. 

For the one loop propagator corrections one has to first compute the corrections to all the 2-point functions which were relevant 
for the derivation of ( 1152b . From [36] this is of the form 

p 2 5 ac _ ff uo d X f xh ac jj f aod X 

-ffb -( p 2 +M 2 )8 ac 8 bd ) \Nf cab QS« c S bd + Wf ace f bde +Rf abe f cde +Sdf cd J K 1 ' ' 

which is written with respect to the basis |-^4^,(p"/' j and as we work in the Landau gauge the common Lorentz structure, 

Pfjv(p), has been factored off. The first matrix corresponds to the tree part of the action and the quantities X, U, N, Q, W, R and 
S represent the one loop corrections and we have used the shorthand coupling constant a = g 2 / (16ft 2 ). The totally symmetric 
object df cd is defined by, JH, 

d fcd = ^ Ti ^ T a T (b T e T d)^ ^ (154) 

where (T£)b c = — if abc is the adjoint representation of the color group generators. At this stage we note that ( 1153b represents a 
formal definition and no vacuum bubble expansion has been performed. To one loop one can formally invert (1153b to obtain the 
one loop corrections to all the propagators (1 152b which is 



(p 2 +M 2 ) c- cp f fcpq 

~f fpcd 1 Scpst/q 1 fcdr fpqr 

' [{p 2 ) 2 +M 2 p 2 +C A f\ J (p 2 +M 2 )° + {p 2 +M 2 )[(p 2 ) 2 +M 2 p 2 +C A f] J 1 

( Ah cp Cf cpq 



\Ef pcd Gb cp b dq +J f cpe f dqe +Kf cde f pqe + Ld c ^ pq ) fl + <^( a ) ■ (\5^) 



The objects A, C, E, G, K, J and L are related to the quantities of the one loop matrix of (11531 ). However, as we are focussing in 
this article on the gluon propagator at zero momentum then we only need the relation for A and note that the formal correction 
at one loop for this is 

A = [(pif+Jpi + crfp x [(p 2 +M 2 ) 2 X-C A7 2 (jV + ^)(p 2 + M 2 )+C A y 4 (6 + C^ + iQW)] . (156) 

As noted above we could in principle compute the exact form of each of the 2-point functions contributing to ( 1156b but 
ultimately as we will take the p 2 — > limit this would be unnecessarily overcomplicated. Instead we compute those pieces of 
( 11561 ) which remain at leading order in the vacuum bubble expansion. 

For this we need to determine the fourteen contributing Feynman diagrams. These were generated using the Qgraf package, 
rf38tl . and converted into Form input language where FORM is a symbolic manipulation language, rf39tl . The vacuum bubble 
expansion written in Form was applied to each integral and expressions obtained for all the 2-point functions. As these depend 
on M 2 and y 2 we were able to check that our expressions agreed with those already determined in the M 2 — 0, y 2 + case 
of lUBt. Moreover, we also checked the explicit Slavnov-Taylor identities for the renormalization of the new mass operator in 
the MS scheme by applying the MINCER algorithm, [40], written in FORM, [41], to the Green's function where the operator 
(ipff'cpjf' — W b uf b ) is inserted in an co 2-point function. The resulting renormalization constants were crucial to not only ensuring 
that our conventions were consistent but also that our 2-point function vacuum bubble expansion is correctly finite after being 
fully renormalized. The upshot of our computations is the observation that for the gluon propagator in the zero momentum limit 
only X is required for the leading (momentum independent) term of (11561 ). Thus we finally obtain 



A, 4 167t 2 A, 8 



M 4 9 /—a m ? M 6 ( 9 A, 4 \ 15 2 , 4 1 3, 4 1 ml 

2 



15, a 1 m 2 , /9 21 A 4 \ 3 /—a m 2 

-n^+M 2 n^r \/M 4 -4A 4 n^ 

3 m\ \% 16 M 4 J 16 v m\ 



8 (VM 4 -4A 4 ) 



(157) 



for the one loop correction at zero momentum where all mass variables correspond to renormalized ones. We note that unlike 
the M 2 = case the nonzero freezing akin to tree order is driven by the gluon 2-point function correction. By contrast in the 
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FIG. 7: The gluon propagator ©O(0) drawn in function of M 2 in units A«s = 1. 



M = situation the gluon suppression at one loop derives from the cp 2-point which is related to the horizon condition and the 
gap equation. Also in this case A will be 0(p 2 ) and not 0(1) to retain suppression at one loop, l36ll . 

We apply the variational principle to the gluon propagator in a completely similar manner as in the case of the ghost propagator. 
Therefore, we replace M 2 with (1 — l 2 )M 2 andg 2 with lg 2 in the expression (I157K expand up to order i 1 , and set I = 1. Doing so, 
we find the original expression (1157b for the gluon propagator back. Firstly, we try to apply the principle of minimal sensitivity. 
Therefore, we have depicted the gluon propagator in FIG. |7] First, we notice that £>W(0) is positive for all M 2 e [0, 1.53 
Unfortunately, we do find neither a minimum nor a point of inflection in this interval. Therefore, we shall take the value of M 
which was obtained in the study of the ghost propagator (see previous section). Hence, setting M 2 = 0.37A^-, gives 



2 



m , . 0.63 11.65 
- (1) (0) = ^ = ^. (158) 

MS 

Evidently, the effective coupling is still smaller than 1, cfr. ( 1151l >. 

In summary, the infrared value of the ghost propagator and the zero momentum gluon propagator seem to be reasonable. We 
find a non-enhanced ghost propagator and a gluon propagator which is non-zero at zero momentum. Our results for the gluon 
and ghost propagator are of a qualitative nature as we are only working in a first order approximation. In order to improve these 
numerical results, higher order calculations are recommendable. This is however far beyond the scope of the present article. 



3. The temporal correlator: violation of positivity 

With the help of the variational technique, we can also show that the gluon propagator displays a violation of positivity. If we 
rewrite the gluon propagator in the Kallen-Lehmann spectral representation, 



lp p 2 +M 2 



^{p 2 ) = f o dM 2 D -^L, (159) 



p(Mp) should be a positive function in order to interpret the fields in terms of stable particles. If p(M 2 ) < for certain M 2 , 
D (p 2 ) is positivity violating. As a practical way to uncover this property, one defines the temporal correlator ||2|] 

C(t) = ' dM p p(Ml)e- M »> = i- J + \-'<*v{p 2 )&p . (160) 

Consequently, if we can show that c(t) becomes negative for certain f, p(M 2 ) cannot be positive for all M 2 , resulting in a 
positivity violating gluon propagator. If the gluon propagator vanishes at zero momentum, © (0) = 0, one can immediately 
verify from ( 1159t that p{M 2 ) cannot be a positive quantity. However, having £>(0) + 0, does not exclude a positivity violation 
as we shall soon find out. 



We can now apply the variational technique on the temporal correlator. At tree level, this c(t) is given by 

J_ r+- c _ ipt P 2 +M 2 
2W-~ p*+M 2 p 2 + {l 2 {M 2 )) 



C(t,M 2 ) = ^/ + V' P ' - ,/' +M ' ? dp, (161) 
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FIG. 8: C(t) (fm) in function of t (fm). 
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FIG. 9: g 2 N/(l(m 2 ) in function of; (fm). 



where X 2 (M 2 ) is still determined by the gap equation ( 11491 >. Replacing M 2 — ► (1 — l 2 )M 2 and g 2 — > £g 2 is redundant in this case, 
as we only have the tree level gluon propagator <D (p 2 ) at our disposal. We shall now implement the minimal sensitivity principle 
as follows: for each different value of f, we minimize the temporal correlator with respect to M 2 . C{t) displays a minimum at 
M^ in * 0, for t > 6/Ass. In TABLE|ffl] some values for M^ in (f ) for different t are presented. For t < 6, we have taken M 2 = 0; it 
is clearly visible from the table below that M^ in — * for decreasing t. The corresponding c(f,M^ in ) is depicted in FIG. [8] Both 





6 7 8 9 10 


M 2 ■ 

mm 


0.16 0.35 0.51 0.65 



TABLE III: Some M 2 - for different t in units A m = 1. 



the x-axis and y-axis are shown in units fm (1/A„s = 0.847 fm), in order to compare our results with j29ll42Tl . Not only do we 
find a positivity violating gluon propagator as c(t) becomes negative, but even the shape of this function is consistent with the 
lattice results 4 ! 29ll42Tl . Moreover, in Il29[|42ll . the positivity violation starts from t ~ 1 .5 fm, in good agreement with our results. 
Finally, FIG. [9] displays the corresponding values of g 2 N/l6% 2 . We can conclude that the previous results are reliable for t < 8 
as g 2 N/16n 2 is smaller than one. 



1 29] included quarks, while |43| considered gluodynamics as we are studying in this work. 
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C. A remark about the strong coupling constant 

A renormalization group invariant definition of an effective strong coupling constant g 2 ^ can be written down from the knowl- 
edge of the gluon and ghost propagators as 

glffip 2 ) = 8 2 (f) V (/>V) Q 2 (P 2 ,f) , (162) 
see e.g. [3]. £> and Q stand for the gluon and ghost form factor, defined by 

i(p 2 ) = p 2 v(p 2 ), 

Q{P 2 ) = P 2 Q{P 2 )- (163) 
The definition ( 11621 i represents a kind of nonperturbative extension of the nonrenormalization of the ghost-gluon vertex. At 
the perturbative level, this is assured by the Ward identity (1541 1. Z g = Z~ l Z A '^ 2 . Usually, this is assumed to remain valid 
at the nonperturbative level. Although this cannot be proven, this hypothesis has been corroborated by lattice studies like ll^Pfl . 

In recent years, there was accumulating evidence that g 2 s (p 2 ) would reach an infrared fixed point different from zero: see 
e.g. JH BIH for a Schwinger-Dyson analysis, ll36i 14511 in the ordinary Gribov-Zwanziger approach and l46l l47tl for lattice 
results. These studies are mostly done in a MOM renormalization scheme, with the exception of ll36tl where the MS scheme 
was employed. The manifestation of this infrared fixed point was motivated in Schwinger-Dyson studies and the ordinary 
Gribov-Zwanziger case by means of the power law behavior of the form factors, 

®Cp 2 W « (p 2 f a > 

£(A*»o ~ (p 2 y a > < i64 > 

being expressable in terms of a single exponent a. The Schwinger-Dyson community heralded in a variety of studies the value 
a ss 0.595, whereas the Gribov-Zwanziger scenario gives a = 1. Anyhow, substituting a behavior like ( 11641 > into the definition 
(1162t leads to g 2 ^{p 2 ) p i~Q x (p 2 ) , opening the door for a finite value. 

However, once again quoting the more recent large volume lattice data of lfT6l[T8l[T9ll43ll . the power law behavior (1164b seem 
to be excluded in favor of 

■D(p 2 ) p 2^ Q oc p 2 , 

i(p 2 W <* (P 2 )° > (165) 

leading to a vanishing infrared effective strong coupling constant at zero momentum since gl s (p 2 ) p 2 a0 « p 2 . The refined 
analysis in this paper of the extended Gribov-Zwanziger action, including an additional dynamical effect, allows us to draw 
a similar conclusion up to the one loop level, i.e. an infrared vanishing g 2 ^. Certain lattice studies also pointed towards this 
particular scenario 114811 . 



V. THE BRST BREAKING IN THE GRIBOV-ZWANZIGER THEORY 

We recall here that the Gribov-Zwanziger action (0 is not invariant under the BRST transformation (122V Indeed, if we take 
the BRST variation of the action (0, one finds a breaking term Ay given by 

A Y = sS = gy 2 J d 4 xf bc (Ay; - (D™c m ) (iff + ) . (166) 

We see that the presence of the Gribov parameter y prevents the action from being invariant under the BRST symmetry. 
Nevertheless, this fact does not prevent the use of the Slavnov-Taylor identity to prove the renormalizability of the theory, 
which is very remarkable. Since the breaking Ay is soft, i.e. it is of dimension two in the fields, it can be neglected in the 
deep ultraviolet, where we recover the usual notion of exact BRST invariance as well as of BRST cohomology for defining the 
physical subspace j49ll . However, in the nonperturbative infrared region, the breaking term cannot be neglected and the BRST 
invariance is lost. In the following, we shall present a detailed analysis of this breaking and of its consequences. In particular, 
we shall be able to prove that the origin of this breaking can be traced back to the properties of the Gribov region Q.. Moreover, 
it turns out that the existence of this breaking enables us to give an elementary algebraic proof of the fact that the Gribov 
parameter y is a physical parameter of the theory, entering thus the expression of the correlation functions of gauge invariant 
operators. 
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A. The transversality of the gluon propagator 

The reader might wonder whether the gluon propagator still remains transverse in the presence of the Gribov horizon. As 
the gluon propagator is the connected two-point function, we ought to consider the generator Z' of connected Green functions, 
which can be constructed from the quantum effective action 5 F by means of a Legendre transformation. The renormalizability 
of the theory entails that T obeys the renormalized version of the Ward identity d32b . or 



Ml • ( 167 ) 



ST 

bb" 

Introducing sources I a {J£) for the fields b a (A^) and performing the Legendre transformation, the identity d 1 67b translates into 

S7 6 

/« = 3^, (168) 
Acting with on this expression, and by setting all sources equal to zero, we retrieve 



8 2 Z 6 

o = d; 



" &/£(*)&*(y) 

which expresses nothing else but the transversality of the gluon propagator. 



= 3J(A^(y)) , (169) 

/,/=o 



B. The BRST breaking and its consequences on the Slavnov-Taylor identity 



Let us present here a few considerations on the consequences stemming from the BRST breaking Ay appearing in the left hand 
side of equation (1166b of the Slavnov-Taylor identity. Our argument will follow 05 Oil . We start from the generalized Slavnov- 
Taylor identity ( TJTb which is fulfilled by the enlarged action E ((28). The quantum effective action, T = E + frT^' + . . ., obeys the 
quantum version of this Slavnov-Taylor identity 1 2811 . 
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(170) 



We now pass to the Gribov-Zwanziger action, defined by giving the sources (M, N, U, V) their physical values d20l and ((24). As 
a consequence, the physical quantum effective action r p h ys will now obey a broken Slavnov-Taylor identity, 
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(171) 



wher eby [ Ay • r p h ys ] represents the generator of the 1PI Green functions with the insertion of the composite operator Ay. Expres- 
sion (1171b generalizes at the quantum level the broken identity of equation ( 1166b . Once having a Slavnov-Taylor identity like 
dm) at our disposal, we can obtain relations between different Green functions by acting on it with test operators g*p( Y| ^ & ¥( Xl! ) » 
with *P any field, and by setting all fields and sources equal to zero at the end. The breaking term in the r.h.s. of expression (117 lb 
will be translated into an extra contribution. In particular, we shall obtain 



5"[5(r phys 



5"[A r r phys ] 



fields, sources=0 



5 v P(jti)...5 v P(>:„) 



(172) 



fields : sources=0 



5 This is the generator of the 1PI Green functions. 
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One sees thus that the r.h.s. of the foregoing expression, corresponding to a 1PI Green function with the insertion of the 
composite operator Ay and with n amputated external legs of the type ^(xi), . . . ^(Xn), gives precisely the modification of the 
relationships among the Green functions due to the Gribov horizon. To our understanding, the contributions stemming from 
the r.h.s. of equation ( 1172b should be correctly taken into account when checking the validity of the Slavnov-Taylor identi- 
ties or when invoking Slavnov-Taylor related identities in computations when the restriction to the Gribov horizon is understood. 

It is worth noticing that the breaking term of (|172t will certainly vanish if the chain ^(xi) . . .^(xn) has a ghost number 
different from +1. Indeed, the action preserves ghost number and the breaking term Ay itself carries a nonvanishing ghost 
charge of +1, so that the operator gip^ giy( r ] must have ghost number —1 in order to allow for a nonvanishing contribution 
( fT72l) . 

In summary, we emphasize that the broken Slavnov-Taylor identity ( 1171t does in fact maintain a powerful predictive character. 
It allows us to establish relationships among various Green functions of the theory in a way which takes into account the presence 
of the Gribov horizon. At the same time, there exist Green functions for which the breaking of the Slavnov-Taylor identity is 
harmless. In particular, this is the case when considering gauge invariant operators built up with only the gauge fields A a ^. For 
these Green functions, the physical quantum action r p t,y S behaves as it fulfills the unbroken Slavnov-Taylor identity, namely 

5(r phys ) = 0. (173) 

The gauge invariance of the correlator (F 2 {x)F 2 (y)) implies in fact that no useful information can be extracted for it from the 
Slavnov-Taylor identities. In a loose way of speaking, (F 2 {x)F 2 (y)) lives on its own and is not related to other Green functions. 
To formally prove this, one should add the operator F 2 (x) to the action with a (BRST invariant) scalar source K(x), the (broken) 
Slavnov-Taylor identity dl7U will remain unchanged. Hence, similarly as in the previous subsection, by performing a Legendre 
transformation to pass to Z' and by acting with on that identity and again setting all sources to zero, it will follow that 

there is a trivially vanishing breaking term due to ghost charge conservation. 



C. A few words on unitarity 

Certainly, the BRST breaking and its consequences on the Green functions of the theory deserve further investigation. In this 
respect one could attempt to evaluate some gauge invariant correlation function like, for instance, (F 2 {x)F 2 (y)) in order to see 
if, despite the presence of the BRST breaking and of a positivity violating gluon propagator, this gauge invariant correlation 
function might displays a real pole in momentum space. A first hint that something like this might happen, has been given by a 
tree level computation in 18J] . 

As one can easily figure out, the presence of the BRST breaking Ay is related to the lack of unitarity in the gluon sector. 
To our understanding, this is a manifestation of gluon confinement: unitarity is jeopardized in the gluon sector because gluons 
are confined. This is also apparent from the positivity violation exhibited by the gluon propagator, which does not allow for 
a physical interpretation of the elementary gluon excitations. One might have the tendency to believe that the existence of the 
soft breaking Ay of the BRST symmetry is a welcome feature, in particular signalling that, in a confining theory, physics in 
the infrared region is not necessarily definable in the same way as in the deep ultraviolet, where the BRST breaking could be 
neglected and one recovers usual perturbation theory. As we already stated in the beginning of this section, in the ultraviolet, we 
also recover the usual notion of the BRST cohomology 04911 . allowing to prove that the ghost degrees of freedom cancel against 
2 unphysical gluon polarizations, leaving over only 2 physical transverse polarizations, endowed with a positive norm. In the 
confining regime, it is unknown what the analogue of this scenario might be. The absence of the BRST symmetry in the infrared 
does not necessarily entail that the theory is not unitary. Certainly, the S -matrix of the excitations of the physical spectrum has to 
be unitary. But as gluons are not the excitations belonging to the physical spectrum, unitarity is not to be expected in the sector 
described by the elementary gluon fields. From this perspective, the question of what the number of physical gluon polarizations 
might be in the nonperturbative confining infrared sector loses its context. 



D. The BRST breaking as a tool to prove that the Gribov parameter is a physical parameter 

The breaking term d!661 l has also the interesting consequence that it allows us to give a simple algebraic proof of the fact 
that the Gribov parameter y is a physical parameter of the theory, and that as such it can enter the explicit expression of gauge 
invariant correlation functions like for instance (F 2 (x)F 2 (y)) or the vacuum condensate (F 2 ). In fact, by taking the derivative 
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of both sides of equation (1166b with respect to y 2 one gets, 

^ = ^=sj dV (a;< - (flJV) (< + tf ) ) , (174) 

from which, keeping in mind that the BRST operator s as defined in equation (f22b is still nilpotent, it immediately follows that 
Jp- cannot be cast in the form of a BRST exact variation, namely 

dS 



^a*sA y , (175) 



for some local integrated dimension two quantity Ay. From equation ( 1175b it becomes then apparent that the Gribov parameter y 2 
is a physical parameter, as much as the gauge coupling constant g, for which a similar equation holds. Furthermore, it is worth 
underlining that, due to the form of the BRST operator s, the presence of the soft breaking Ay is, in practice, the unique way to 
ensure that the Gribov parameter indeed is a physical parameter and not an unphysical one, as it would be the case of a gauge 
parameter entering the gauge fixing term. Let us suppose that the part of the action Sy containing the Gribov parameter would 
be left invariant by the BRST transformation ( 1221 ). namely 

sS y = 0, (176) 

instead of inducing the breaking term Ay. Since Sy depends on the auxiliary fields fef, <p", c , C0^ e , co^ which constitute a set of 
BRST doublets 6 $2^, it would follow from equation ( 11761 ) that a local integrated polynomial Sy would exist such that 

Sy = sSy. (177) 

Subsequently, taking the derivative of both sides of expression ( 1177b with respect to y 2 , one would obtain 

dSy dSy 

3y 2 ' 3y 2 ' 

a relation implying that y 2 would have the same meaning as an unphysical gauge parameter 7 . In turn, this would imply that 
correlation functions of gauge invariant operators would be completely independent from y 2 . We see thus that the presence of 
the soft breaking term Ay plays an important role, ensuring that y 2 is a relevant parameter of the theory. The same conclusion 
also holds when the Gribov-Zwanziger action is supplemented by the BRST invariant mass term doTT i. The existence of the 
breaking Ay thus seems to be an important ingredient to introduce a nonperturbative mass gap in a local and renormalizable way. 

A question which arises almost naturally is whether it might be possible to modify the BRST operator, i.e. s — > s m , in such a 
way that the new operator s m would be still nilpotent, while defining an exact symmetry of the action, s m S" = 0. Although we are 
not going to give a formal proof, we can present a simple argument discarding such a possibility. We have already observed that 
the BRST transformation ( 1221 defines an exact symmetry of the action when y = 0, which corresponds to the physical situation 
in which the restriction to the Gribov region has not been implemented. Hence, it appears that one should search for possible 
modifications of the BRST operator which depends on y, namely 

s,„ = s + sy, (179) 

whereby 

Sy = y-dependent terms , (180) 

so as to guarantee a smooth limit when y is set to zero. However, taking into account the fact that y has mass dimension 
one, that all auxiliary fields ^p^ e ,9^ c ,G^f, CO™'^ have dimension one too, and that the BRST operator s does not alter the 
dimension of the fields 8 , it does not seem possible to introduce extra y-dependent terms in the BRST transformation of the fields 
( (p^ c , (p^ c , CO™ , C0^' J while preserving locality, Lorentz covariance as well as color group structure. 



6 We remind here that a BRST doublet is given by a pair (a, P) transforming as: .sa = p, sp =0. It can be shown that a BRST doublet has always vanishing 
cohomology, meaning that any invariant quantity, sF(a, p) = 0, has necessarily the form of an exact BRST cocycle, namely F(a, P) = sF(a, p). 

7 One easily shows that in this case, = — {^^jp-§^=0 for any gauge invariant operator § . 

8 It is understood that the usual canonical dimensions are assigned to the fields A",, b", c" , c" 1231 . as shown in TABLEU It is apparent that the BRST operator 
* does not alter the dimension of the fields. 
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E. Tracing the origin of the BRST breaking 

Having clearly seen the explicit loss of the BRST symmetry, it would be instructive to point out more precisely where this 
breaking originates from. We recall that the BRST transformation of the gluon field A p is in fact constructed from the in- 
finitesimal gauge transformations. Indeed, for an infinitesimal gauge parameter co", the corresponding gauge transformation is 
determined by 

8^ = Df<d b , (181) 

which can be compared with the BRST transformation (|22| |. Based on this identification, we shall present our argument using 
infinitesimal gauge transformations. In particular, we shall establish the following proposition: any infinitesimal gauge transfor- 
mation of field configurations belonging to the Gribov region CI, necessarily gives rise to configurations which lie outside of £2. 
We can distinguish 2 cases. 

• The field is not located close to the boundary d£l 

Let us consider a gauge configuration A M which belongs to the Gribov region £2 but not close to its boundary dQ. (the 

horizon), thus d M A M = and — 3 /U D 1 ,(A) > 0. Next, consider the field A p obtained from A^ through an infinitesimal gauge 
transformation with parameter co, 

A„=A tl +D ll (A)<B. (182) 

This configuration A^ cannot belong to £2. Suppose the contrary, then d^A^ = = d^A^ would lead to 

8 A ,D„(A)co = 0, (183) 

in contradiction with the hypothesis that A^ is not located on the boundary 3£2, thus there are no zero modes co allowing 
for (fT83l to hold. 

• The field A^ is located close to the boundary 3£2 

In this case, we can even make a more precise statement. If A A , lies very close to the boundary 3£2, we can decompose it as 

A,i = a^ + C^, (184) 

with Cf, G d£l, thus C p lies on the horizon. The shift is a small (infinitesimal) perturbation. Obviously, d^C^ = d^a^ = 0. 
Subsequently, we find 

Ap, = C p + a/J + D M (C)(i)+... (185) 

for the gauge transformed field at lowest order in the infinitesimal quantities CO and a^. Since Cu E 3£2 and by identifying 
co with the zero mode corresponding to C^, we find 

M,, = dMC)(0 = . (186) 

showing that A^ is transverse. The field A^ also lies very close to the boundary 8£2. However, as it follows from Gribov's 
original statement 9 10], it is located on the side of the horizon opposite to that of the field A^, i.e. it lies outside of the 
Gribov region £2. 

We can conclude thus that any infinitesimal transformation of a gauge field configuration which belongs to the Gribov region 
fi, results in another configuration which lies outside Q.. Since the BRST transformation of the gluon field is naturally obtained 
from the infinitesimal gauge transformations, it is apparent that the breaking of the BRST symmetry looks almost as a natural 
reflection of the previous result. 

We can also offer a pictorial depiction of what is happening. We recall that the Gribov region £2 is convex, bounded in all 
directions in field space, that every gauge field has an equivalent representant within £1, that the origin A^ = belongs 10 to 



9 For the benefit of the reader we quote here Gribov's statement, proven in | 7]: for each field A,, belonging to the Gribov region Q. and located near the boundary 
dQ., i.e. A fI = C l j + a^, there exists an equivalent field A I1 ,A I1 = +a„+D„(C) 09, near the boundary 3£1, located, however, on the other side of the horizon, 
outside of the region SI. 

This means that perturbation theory belongs to SI. 
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n SHI H3], and that every gauge configuration near the horizon 8£2 has a copy on the other side of dQ. 0]. The first 4 
quoted properties are important to make Q. a suitable domain of integration in the path integral, i.e. we can restrict the whole 
space of A A ,-configurations to Q. as proposed by Gribov. However, implementing this restriction in A^-space jeopardizes the 
BRST invariance. As we have seen, if we move throughout A^ space with a BRST transformation (cfr infinitesimal gauge 
transformations), we must unavoidably cross the horizon dd. Hence, restricting the fields within the horizon breaks the BRST 
invariance. 



F. The Maggiore-Schaden construction revisited 

The authors of the paper [23] attempted to interpret the BRST breaking as a kind of spontaneous symmetry breaking. We 
shall now re-examine this proposal and conclude that, instead, the BRST breaking has to be considered as an explicit symmetry 
breaking, where we shall present a few arguments which have not been considered in 12311 . Although this discussion might 
seem to be only of a rather academic interest, there is nevertheless a big difference between a spontaneously or explicitly broken 
continuous symmetry, since only in the former case a Goldstone mode would emerge. For the benefit of the reader, we shall first 
explain in detail the approach of ll23"ll . One starts by adding the following BRST exact term to the Yang-Mills action: 

Si = sf d 4 x(lfd^* + (3?dvD?<? b ;) , (187) 

with s, the same nilpotent BRST operator as defined in d22l >. The first term represents the Landau gauge fixing, while the second 
term is a BRST exact piece in the fields ((p,co,(p,Co). Of course, from expression ( |187t , it follows that s defines a symmetry 
of the action Sym + S\. As a consequence, the nilpotent operator s allows us to defines two doublets (<p,0>) and (9,©). This 
doublet structure implies that we can exclude these fields from the physical subspace j28il49ll . which makes Sym +Si equivalent 
to the ordinary Yang Mills gauge theory. Next, Maggiore and Schaden introduced a set of shifted fields, which -translated to our 
conventions- are given by: 

< = <" + Y 2 5 £ %, 



C + KT T CO,, X. 



•ft j 



b a = l/ a +g>ff abe tffx M . (188) 
All fields ((plf 6 ,<p!f 6 ,c' a ,b' a ) have vanishing vacuum expectation value (VEV), namely 

W?) = tf a ) = (b' a )=0. (189) 
Along with these new fields (q#* 9jf , c'\ b la \ one introduces a modified nilpotent BRST operator s given by: 

sc' cl = b' a , Jb' a = , 

IA C ; = -Dfc b , 3fflJ* = 0, (190) 

which looks exactly like ( l22b . However, we emphasize that by introducing these new fields, the BRST operator Fwill give rise 
to an explicit x-dependence when acting on the field cd^f: 



s(0 



■fl Tfi 



?J? 6 +Y 1 S a V (191) 



Furthermore, by taking the vacuum expectation value of both sides of equation dl9U . one gets 

(7<)=y 2 5 fl %, (192) 



from which the authors of M23I1 infer that the BRST operator s suffers from spontaneous symmetry breaking. Notice also that 
( |192t gives a VEV to a quantity with a free Lorentz index. 

With the introduction of the shifted fields, we can rewrite the action Si as: 

Si = r|d 4 Jc(c^9^ + o52 c 9vZ)f^^ . (193) 
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The last two terms can be simplified, leading to 

5! = ijd'x^d^ + Wp.Dt^-g^fU^ 



(194) 



If we calculate this action explicitly, we recover the original Gribov-Zwanziger action, without the constant part 4-^(N 2 — 1). 
For this reason one adds — y 2 JJ d 4 Jc3 v co^" to the action Si. Doing so, one finds 



S\ = s 



= /d 4 x 



b' a d^A;+c' a d 



I 



+ d*X 



d 4 x 



(195) 



-gf^ffabc^ - g^'fatcD^c* - 4y> 2 - 1 ) 
If we naively assume that we can perform a partial integration, we find after dropping the surface terms 

Sym + Sx = ©- grV^I d 4 x«c rf . 



(196) 



The last expression reveals that one has recovered the Gribov-Zwanziger action from an exact F-variation with the addition of 

an extra term ( —gf"f abc J d 4 xWff'D'ff(r \ . However, this term is irrelevant as we shall explain now. Assume that we want to 

compose an arbitrary Feynman diagram without any external co leg and thereby using the action ( 11961 ). The second term from 
this action can never contribute to this Feynman diagram as it contains an external co. Indeed, this leg requires an co-leg, which 
in its turn is always accompanied by an CO leg. Hence, the action J 196b is equivalent to the standard Gribov-Zwanziger action ^ 
when we exclude the diagrams containing external co legs 11 . 

Although at first sight this construction might seem useful, it turns out that a few points have been overlooked. Let us 
investigate this in more detail. Firstly, we point out that rather delicate assumptions have been made concerning the partial 
integration. To reveal the obstacle, we perform once more the partial integration explicitly, 



Jd 4 xfx^ v Df cpf ' = surface term - J d^xfS^Df^ 



Iba 



(197) 



Normally, one drops the surface terms, as the fields vanish at infinity. However in this case, as Xf, does not vanish at infinity, it is 
not sure if the surface terms °< x p will be zero. One would have to impose extra conditions on the fields to justify the dropping 
of the surface terms. On the other hand, when we do not perform the partial integration to avoid the surface terms, we are facing 
an explicit, unwanted x-dependence in the action, resulting in an explicit breaking of translation invariance. 
Another way of looking at the problem consists of performing a partial integration on the second term of the action d 1 94-t before 
applying the BRST variation J. Doing so, we find, 



S\ = s 



d 4 x (e%A£ - avco^fcpf ' - gfafifabcAl - y 2 ^) 



(198) 



Subsequently, applying the BRST variation gives, 



Si 



d 4 x 



d 4 * 



d 4 * 



© - 81 2 f' bc f ^xafD^'c 11 . 



(199) 



In this case, we do not encounter the problem of nonvanishing surface terms. To recapitulate, if we first let the BRST variation 
act on the action ( 1194b . and then perform a partial integration, we find a different result than performing these two operations 
the other way around. This difference is exactly given by the surface term from equation ( 11971 ). This discrepancy arises of 
course from the explicit x-dependence introduced in the BRST transformation J, giving nontrivial contributions. For example, 



For our purposes these diagrams are irrelevant, e.g. the vacuum energy, the gluon and ghost propagator, 
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we introduced a term — y 2 ?J d 4 x3 v co"" which might seem to be zero since we are looking at the integral of a complete derivative 
(thus usually taken to be a vanishing surface term), but when the BRST variation is taken first, a nontrivial integrated piece 
remains. 

Apparently, to find the correct Gribov-Zwanziger action with the Maggiore-Schaden argument, there is some kind of a "hidden 
working hypothesis" that J 198b is the correct action to start with, and that partial integration is not always allowed 12 . The 
fact that there seems to be a kind of "preferred" action to start with, is just a signal that there is a problem with the boundary 
conditions for some of the fields and hence surface terms when integrating. 

Even if one forgets about the previous criticism, a second problem arises. In the Gribov-Zwanziger approach, we recall that 
the parameter y is not free and is determined by the horizon condition ( TToT l. As it has been explained in section UTaI the solution 
y = is excluded. In equation ( 11351 *, we found a solution for y + 0. This gave rise to a positive vacuum energy E vac > 0, as one 
can see from equation (11371 >. see also 12511 . However, according to Maggiore-Schaden argument, at one loop order the stable 
solution should be that corresponding to y = ll23ll . as, if y = 0, the vacuum energy would be vanishing, i.e. E vac = 0, which is 
energetically favored over a positive vacuum energy. This delivers a contradiction with the Gribov-Zwanziger approach, as the 
restriction to the Gribov region requires that y + 0, thus giving a positive energy E vac > at one loop. 

To end this section, let us now consider the new operator / d 4 x (yfqf - off <df} within the Mag giore-Schaden approach. 
We observe that we obtain an explicit x-dependence if we rewrite this operator in terms of the new fields, 

J d 4 x (cpf cpf - «) = J d 4 x (q>';% ah - « + ^x^'"" - y^r - j\^(N 2 - 1 )) . (200) 

However, this x-dependence is necessary so that (I2001 > would be invariant under the new BRST symmetry I, 

I J d 4 x (cp^cp^ - « + Y%<fC - Y^-VPr - l\x,{N 2 - 1 )) = . (201) 

A second option is to introduce the BRST F-exact mass operator 

* J d 4 x (iEgV) = J d 4 x (y'fyf - cpf cpf + fx^f) , (202) 

which also displays an explicit x-dependence. 

Finally, let us consider a third and last possible option. If we would have started with the following mass operator, 

|d 4 x(^< & -«) , (203) 

which does not contain an x-dependence, this operator is not left invariant by the symmetry J. In fact, 

sJd 4 x(y'£ b tf b -Wft<of^ = -f j d A x Xfl <tf? ■ (204) 

One sees that with the introduction of the new mass operator, the Maggiore-Schaden construction will always give rise to an 
explicit breaking of translation invariance if the BRST invariance J has to be preserved. We thus conclude that the Maggiore- 
Schaden construction cannot be implemented in the presence of the new operator and even without the new mass operator we 
have collected a few arguments from which the frame of a possible spontaneous symmetry breaking cannot be applied to the 
Gribov-Zwanziger action. 



G. A few remarks on the Kugo-Ojima confinement criterion 



In this section we shall take a closer look at the Kugo-Ojima confinement criterion M24H in relation to the Gribov-Zwanziger 
action. In the literature, it is usually stated that the Kugo-Ojima confinement criterion is realized when the Gribov-Zwanziger 
scenario is realized. A key ingredient in the criterion is m(0) = —1, whereby u(Q) is the value at zero momentum of a specific 
Green function, u is related to the ghost propagator in the Landau gauge according to 115 ill 

g(p 2 ) p 2^ = \— ^r. (205) 



If it would be allowed, one would be able to cross from the second action 11981 to the first one U94L but as we have just shown, these two starting actions are 
inequivalent. 
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From this expression, it is obvious that an infrared enhanced ghost propagator results in m(0) = — 1, thereby fulfilling the criterion. 
Let us recall here that the derivation of the Kugo-Ojima criterion is based on the assumption of an exact BRST invariance and is 
written down in a Minkowskian rather than an Euclidean space-time. This has a few repercussions: 

• At a nonperturbative level, some care should be taken when passing from Euclidean to Minkowski space-time. According 
to our understanding, it is not clear whether a Wick rotation can always be implemented. E.g., the gluon propagator (1 1 031 > 
can exhibit two complex conjugate poles, so one should be careful of not crossing these poles when the contour is Wick 
rotated. Clearly, there could be potential caveats when considering a more complicated gluon propagator. 

• A more crucial shortcoming is the following. As we have emphasized in the foregoing section, the restriction to the Gribov 
region inevitably leads to a breaking of the BRST symmetry which, however, was the very starting point of the Kugo- 
Ojima analysis. In addition, parts of the Kugo-Ojima study rely on analyzing the charge of the global color current and 
the expression of a piece of it in terms of the BRST symmetry generator. In our opinion, as the Gribov-Zwanziger action 
is essentially different from the usual Faddeev-Popov fixed (Landau gauge) action due to new fields, extra interactions 
and especially another symmetry content, the Kugo-Ojima analysis cannot simply be applied to the Gribov-Zwanziger 
formalism, although both might superficially seem to be in accordance with each other. Therefore, it seems to us that one 
cannot verify the Kugo-Ojima criterion (u(0) = — 1) when the restriction to the Gribov horizon is taken into account 13 . 

• The latest lattice data point towards a ghost propagator which is no longer enhanced, so that the condition u (0) = — 1 does 
not seem to be realized anyhow. 



VI. DISCUSSION 



Our starting point was the original localized Gribov-Zwanziger action, Sqz, and the observation of the new lattice data, which 
shows an infrared suppressed, positivity violating gluon propagator, nonvanishing at the origin and a ghost propagator which 
is no longer enhanced. However, the propagators corresponding to the original Gribov-Zwanziger action are not in accordance 
with these new lattice data. Hence, we have searched for a solution by looking at nonperturbative effects like condensates. 

Therefore, we have added two extra terms to the Gribov-Zwanziger action, Sm = M 2 J d 4 x (<p<p — coco) + 2 ^ A/ 1 - > qk 2 . A first 



intuitive argument why we added the first term, M 2 J d 4 x ((pep — coco), was the following. In the Gribov-Zwanziger action, Sgz, 
an Aip-coupling is already present at the quadratic level. Therefore, altering the cp-sector will be translated to the A-sector, thus 
modifying the gluon propagator. Secondly, this condensate is already present perturbatively as, at lowest order, we have found 

((pep -coco) = 3 ^ ~^ A, 2 , (206) 

with A- 4 = 2g 2 Ny 4 . This implies that the condensate is nonvanishing for y + already in the original Gribov-Zwanziger action. 
It was therefore very natural to add this operator to the theory. The second pure vacuum term, M 2 J d 4 x 2 ^ lN l \ k 2 , was added 
in order to stay within the horizon or equivalently, to keep o(0) smaller than 1 when the horizon condition is implemented. We 
have fixed q by imposing ^ =0; this ensures a smooth limit to the original Gribov-Zwanziger action. 

The extended Gribov-Zwanziger action, Sgz + Sm, has many interesting features. Not only is this action renormalizable, it 
is also remarkable that no new renormalization factors are necessary for the proof of its renormalizability, meaning that only 
two independent renormalization factors are required. As an extra feature, we have also shown that Sgz + Sm + S A 2, with 
S A i = tjt / d 4 xA 2 , is renormalizable. 

Another important observation is that the gluon propagator is already modified at tree level. We have found 

i 2 +M 2 
p 4 + (M 2 +m 2 )p 



This type of propagator is in qualitative agreement with the most recent lattice data, which was the starting point of our analysis. 
In section IV the gluon propagator at zero momentum was also presented at one loop (see equation ( 11571 )). where we switched 



This would also include Schwinger-Dyson results which implemented the restriction to the Gribov region by suitable boundary conditions. 
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off the effects related to A 2 by setting m 2 = 0. By virtue of the novel mass M 2 , £> (p 2 ) + at zero momentum. Also the one loop 
ghost propagator is modified. At small momenta we have obtained, 

$(P%~0 = J2T~0 ' (208) 

with 

<""> - l «fwK* lt * i:ii )- h (" , - , HI-(S)F < 209 » 

We see that the ghost propagator is clearly no longer enhanced, again in accordance with the most recent lattice data. 

Up to this point, the mass M 2 was put in by hand. However, we have treated ((f^'cp" — K>"co") as a composite operator coupled 
to the source J — M 2 . In this way, we have been able to find nonperturbative effects induced by this composite operator without 
altering the original Gribov-Zwanziger action, and making the mass M 2 dynamical. We have developed two methods to find such 
nonperturbative effects. The first method uses the well known principles of the effective action formalism. Unfortunately, the 
calculations become intractable. Therefore, we have implemented a second method, the variational principle. Intuitively, we have 
included effects of the mass term without altering the original Gribov-Zwanziger action by performing a suitable resummation. 
With the help of this technique, we have found in the MS-scheme that a(0), the one loop correction to {j? 2 Q (p 2 )) 2~ is given 
by 14 

o(0) = 0.93, (210) 

resulting in a non-enhanced ghost propagator. Simultaneously, for the one loop gluon propagator at zero momentum, we have 
found 

A|f GeV 2 

which is nonzero. The corresponding value for the coupling constant is smaller than 1, see equation ( 11511 ), which is acceptable 
for a perturbative expansion. We have also checked the positivity violation of the gluon propagator with the help of the variational 
technique and again, our results were in nice agreement with lattice results: not only is the shape of the temporal correlator C(t), 
displayed in FIG.[8]in qualitative agreement, also the value of the point, t ~ 1.5 fm, at which the violation of positivity starts is 
consistent with the results reported in lattice investigations. Using the plots displayed in lfl7ll which were also obtained in the 
SU (3) case, one can extract a rough lattice estimate for the quantities (121 lb and ( 1210b . 

^lattice/Q) 3 (212) 

GeV 2 



P 2 g(p 2 )^0 ~ 5 & o(0) lattice ~ 0.8, (213) 

We notice that our lowest order approximations (121 It and (12101 ) are qualitatively compatible with the current lattice values. 

To conclude, we would like to emphasize that the original Gribov-Zwanziger action already breaks the BRST symmetry. Due 
to this breaking, it is in unclear at present how to define the observables of the theory in the nonperturbative infrared region. 
According to our understanding, this breaking cannot be interpreted as a spontaneous breaking, according to the proposal of 
rf23"tl . In fact, we have argued that the BRST breaking is a natural consequence of introducing the restriction to the Gribov 
region. In addition, we have underlined that the presence of the BRST breaking term in the Gribov-Zwanziger action provides a 
consistent way to ensure that the restriction to the Gribov region can have physical consequences, i.e. that the Gribov parameter 
y enters the expectation value of physical, gauge invariant correlators. In the absence of such a breaking term, the Gribov mass 
parameter would play the role of an unphysical gauge parameter. The presence of the breaking is thus a necessary tool within the 
Gribov-Zwanziger approach, allowing for the introduction of a nonperturbative mass parameter in a local and renormalizable 
way. Finally, we have also commented on the Kugo-Ojima confinement criterion. Since it is fundamentally based on the concept 



We set Ahj- = 0.233GeV, the value reported in |53l . 
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of an exact BRST symmetry, it cannot be straightforwardly related to the Gribov-Zwanziger framework due to the breaking. 

In summary, this paper presented the 4D analysis of the gluon and the ghost propagator within the Gribov-Zwanziger frame- 
work. By comparing these results with recent lattice data, we have found a good qualitative agreement. The ghost and gluon 
propagator have also been extensively studied on the lattice in 2 and 3 dimensions !116lll8Ul9ll53ll . The 3D and 2D analysis of 
the extended Gribov-Zwanziger action, and a comparison with the lattice data, is currently under consideration. 
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APPENDIX: THE ONE LOOP EFFECTIVE POTENTIAL 



We explain in detail how we obtained equation ( 11091 ). We start by evaluating the integral appearing in W(J). We recall that 
this integral originates from: 



dA^ex P -i/d- l u;;(A;;;iL\;; 



det 



a 



-1/2 



= e -4Trln2$ 



(A.l) 



with Q^y = — (d 2 + gf ^ 8^ — dfjd v (l — ~). From this expression it follows that we need to calculate ^Tr lng— to obtain 
the energy functional: 



N z -1 



■ 2g 2 Nf \ 
d 2 -M 2 J 

(d-l) {Tr In (-3 2 (3 2 - M 2 ) - 2g 2 Ny 4 } - Tr ln(-3 2 + M 2 ) } . 



-Trine- = _ (< /-l)Trta(-*-gf-^j 

N 2 -l 



The second part is a standard integral and evaluated as: 

Trln(~a 2 + M 2 



-r(-<*/2) 1 



(4n) d / 2 (M 2 )- d / 2 ' 

with r the Euler Gamma-function. Using dimensional regularization, d = 4 — e we obtain, 



JV 2 -1 



'-l)Trln(-3 2 +M 2 ) = -3 



jV 2 - 1 



5 2 



M 



6 e 



2 v " '' v - 64jt 2 

We recall that we work in the MS scheme. Next, we try to convert the first part in to the standard form, 

N 2 -l 



(d- l)Tr ln(— 3 2 (3 2 — M 2 ) - 2g 2 Ny 4 ) 



N z -1 



N z -1 



1 )Tr In (-3 2 + mfj + Tr In (-3 2 + m\) 



'-1) 
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where we have used the notational shorthand ( 1108b . We still have to calculate the first and the second term of ( 1107b . For the first 
term, we recall that 

A 9 A 9 3 e 2 N 1 

J 4 =Z 2 2 y 4 , with Z 2 2 = l + -i^-, (A.6) 

with Zy2 defined in (55[ . so we find 

~d(N 2 ~l)y 4 = -4(N 2 ~l)y A -4 3 -(N 2 -l)^ l -y 4 + 3 -^y\N 2 -l). (A.7) 
The second term is invariant under renormalization and therefore given by 

From equation dA.41 >. (IA.5t and iA.71 we see that the infinities cancel out nicely, so that the functional energy reads, 

/ \ \ / 4(N 2 -1), 4 d(N 2 -i) 3(JV 2 -l)/8, 4 . m 2 , ml , J\ 

W*\J) = -^^XV^^aV+^i^ + ^ln^l+^ln^-^ln^j , (A.9) 

which is exactly expression ( 11091) . 
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